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Abstract 

We show that in the presence of massive particles such as nucle- 
ons, the standard low energy expansion in powers of meson momenta 
and light quark masses in general only converges in part of the low 
energy region. The expansion of the scalar form factor <r(t), for in- 
stance, breaks down in the vicinity of t = 4M^. In the language of 
heavy baryon chiral perturbation theory, the proper behaviour in the 
threshold region only results if the multiple internal line insertions 
generated by relativistic kinematics are summed up to all orders. We 
propose a method that yields a coherent representation throughout 
the low energy region while keeping Lorentz and chiral invariance ex- 
plicit at all stages. The method is illustrated with a calculation of the 
nucleon mass and of the scalar form factor to order p 4 . 
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1 Introduction 



The effective low energy theory of the strong interaction is based on a si- 
multaneous expansion of the Green functions of QCD in powers of the exter- 
nal momenta and of the quark masses ("chiral expansion"). In the vacuum 
sector, where the only low energy singularities are those generated by the 
Goldstone bosons, dimensional regularization yields homogeneous functions 
of the momenta and Goldstone masses, so that each graph has an unam- 
biguous order in the chiral expansion. In the sector with baryon number 1, 
however, the low energy structure is more complicated. The corresponding 
effective theory can be formulated in manifestly Lorentz invariant form ||T|, 
but it is not a trivial matter to keep track of the chiral order of graphs con- 
taining loops within that framework: The chiral expansion of the loop graphs 
in general starts at the same order as the corresponding tree graphs, so that 
the renormalization of the divergences requires a tuning also of those effec- 
tive coupling constants that occur at lower order - in particular, the nucleon 
mass requires renormalization at every order of the series. 

Most of the recent calculations avoid this complication with a method 
referred to as heavy baryon chiral perturbation theory (HBXPT) [|, || f|, |5| . 
The starting point of that method is the same Lorentz invariant effective La- 
grangian that occurs in the relativistic approach. The loop graphs, however, 
are evaluated differently: The Dirac-spinor that describes the nucleon degrees 
of freedom is reduced to a two-component field and the baryon kinematics 
is expanded around the nonrelativistic limit. At the end of the calculation, 
the amplitude may then be recast into Lorentz invariant form. This method 
keeps track of the chiral power counting at every step of the calculation, at 
the price of manifest Lorentz invariance. 

Mojzis and Fettes et al. |7] have evaluated the irN scattering amplitude 
to order p 3 within that framework. The explicit result is remarkably simple, 
the one loop graphs being expressible in terms of elementary functions. On 
general grounds, the outcome of this calculation must be the same as what 
is obtained if the representation of the scattering amplitude given in |IJ is 
expanded around the nonrelativistic limit - that representation also holds to 
order p 3 . Indeed, Ellis and Tang M have shown that this is the case. 

Quite apart from the fact that the nonrelativistic expansion turns the ef- 
fective Lagrangian into a rather voluminous object and that care is required 
to properly analyze the corresponding loop graphs^, HBXPT suffers from a 
deficiency: The corresponding perturbation series fails to converge in part 
of the low energy region. The problem is generated by a set of higher order 

^^Wave function renormalization, for instance, is not a trivial matter [0, O, O, O]. 
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graphs involving insertions in nucleon lines. A similar phenomenon also ap- 
pears in the effective field theory of the iVW-interaction ||13|| . It arises from 
the nonrelativistic expansion and does not occur in the relativistic formula- 
tion of the effective theory. 

The purpose of the present paper is to present a method that exploits 
the advantages of the two techniques and avoids their disadvantages. In sec- 
tions 3 and 4, we demonstrate the need to sum up certain graphs of HBXPT, 
using the example of the scalar nucleon form factor. Next, we show that 
the infrared singularities of the various one loop graphs occurring in the chi- 
ral perturbation series can be extracted in a relativistically invariant fashion 
(sections 5 and 6). The method we are using here follows the approach of 
Tang and Ellis |J. There is a slight difference, insofar as we do not rely on the 
chiral expansion of the loop integrals - that expansion does not always con- 
verge. A comparison with HBXPT is given in section 7, where we also show 
that our procedure may be viewed as a novel method of regularization, which 
we call "infrared regularization" . The method represents a variant of dimen- 
sional regularization. While that scheme permits a straightforward counting 
of the powers of momentum at high energies, infrared regularization preserves 
the low energy power counting rules that underly chiral perturbation theory. 
Renormalization is discussed in section 8. We then analyze the convergence 
of the chiral expansion of the one loop integrals and show that the expansion 
coefficients relevant for tcN scattering can be expressed in terms of elemen- 
tary functions (section 9). Whenever that expansion converges, the result 
agrees with the one obtained within HBXPT. We also construct an explicit 
low energy representation of the triangle graph, for which the nonrelativistic 
expansion fails. The method is illustrated with an evaluation of the chiral 
perturbation series for the nucleon mass, the wave function renormalization 
constant and the scalar form factor. The physics of the result obtained for 
the form factor and for the cx-term is discussed in sections 12 and 13. We 
demonstrate that a significant part of those infrared singularities that are 
proportional to is common to the form factor and to the irN scattering 
amplitude and thus drops out when considering the predictions of the theory. 
The observation leads to a specific reordering of the chiral perturbation series 
that reduces the matrix elements of the perturbation quite substantially. The 
effects generated by the A(1232) are discussed in some detail (appendix [D]) 
and we also compare our framework with the ancestor of the effective theory 
described here: the static model. Section 14 contains our conclusions. 
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2 Effective Lagrangian 



The variables of the effective theory are the meson field U(x) G SU(2) and 
the Dirac spinor ip(x) describing the degrees of freedom of the nucleon. The 
effective Lagrangian contains two pieces, 

The first part is the well-known meson Lagrangian, which only involves the 
field U(x) and an even number of derivatives thereof. For the second part, 
which is bilinear in ip(x) and if)(x), the derivative expansion contains odd as 
well as even terms: 

r - H 1 ) i r( 2 ) _i_ r( 3 ) _l 

The explicit expressions involve the quantities u, u^, and \±- m the 
absence of external vector and axial fields, these are given by 

u 2 = U , u ti = iu ] dJJu ] , r M = \[u\ d^u] , x± = ui X ui ± u X j u . 

At leading order, the effective Lagrangian is fully determined by the nucleon 
mass and by the matrix element of the axial charge (m and g& denote the 
corresponding leading order values and = + T^): 

£ ( n = $ (ip- + \9A^il^ , 

We disregard isospin breaking effects. The Lagrangian of order p 2 then con- 
tains four independent coupling constants^ 

= Cl ( X+ ) i> i> - £^ Ku u ) $ d»d v iJ> + h.c.) (i) 
+ y i> i> - j i> 7 M i v K, uu] i> . 

Below, we apply the machinery to the scalar form factor. This quantity does 
not receive any contribution from C$ , but two of the terms in C$, 

4 ) = -^(x + )¥^ + |(x + )n(^^) + ... , (2) 

do generate contributions proportional to t\M^ and e^M 2 1, respectively. 

2 We use the conventions of ref. ||. In this notation, the coupling constants of ref. |jj are 
given by: mc° ss = F 2 c\, mc2 SS — —F 2 C4, mc^ ss = —2F 2 cs, m{cf ss — 2 m cf ss ) — — F 2 C2 
(to order p 2 , the terms cf ss and c§ ss enter the observables only in this combination) , while 
those of ref. |Q read: 16 a\ = 8 m c^,+g 2 A , 8 a 2 = 4 m ci— g\, aa = mci, Aa$ = &mci+1— g\. 
In the numerical analysis, we work with F^ = 92.4 MeV, g A — 1-267, M n — M n +. 
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3 Scalar form factor 



We first wish to show that, in the sector with baryon number 1, the standard 
chiral expansion in powers of meson momenta and quark masses converges 
only in part of the low energy region. For definiteness, we consider the scalar 
form factor of the nucleon in the isospin limit (m u = rrid = m), 

(N(P',s')\m(uu + dd)\N(P,s)) =u'ua{t), t = (P' - P) 2 . 

The first two terms occurring in the low energy expansion of this form factor 
were worked out long ago, on the basis of a one loop calculation within the 
Lorentz invariant formulation of the effective theory Jj]]. In that expansion, 
t, rh and M 2 are treated as small quantities of 0(p 2 ), while the nucleon mass 
represents a term of 0(p°). In view of the quark mass factor occurring in 
the definition of cr(t), the low energy expansion starts at order p 2 , with a 

(2) 

momentum independent term generated by C N . 

a(t) = -4 C1 Ml + MMlpL S {t _ 2 Ml) 7 (t) - -^U + 0(p*) (3) 

The constant c\ occurring here is a renormalized version of the bare coupling 
constant in eq. ([!]). Since the renormalization depends on the framework 
used, we do not discuss it at this preliminary stage. The contribution of order 
p 3 is generated by the triangle graph shown in fig. 1 and is fully determined 
by F n and qa- The term involves the convergent scalar loop integral 




Figure 1: Triangle graph. The solid, dashed and wiggly lines represent nu- 
cleons, pions and an external scalar source, respectively. 

1 f d 4 k 1 

7W (2vr) 4 (M 2 — k 2 -ie) (M 2 -(k-qf-it) (m 2 - (P-k) 2 -it) ( ' 

Here and in the following, we identify the masses occurring in the loop inte- 
grals with their leading order values, M w — > M, m N — > m. 
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The function represents a quantity of 0(l/p). Since the external 
nucleon lines are on the mass shell, it exclusively depends on t = q 2 , M and 
m. The function is analytic in t except for a cut along the positive real axis, 
starting at t = AM 2 . The triangle graph also shows up in the analysis of the 
ttN scattering amplitude to one loop, so that the function j(t) is relevant 
also for that case. 

The imaginary part of ^(t) can be expressed in terms of elementary func- 
tions 0: 

Im7(t) = »(' - > ,, ta „ ^gg3 . (5) 

8ny/t (Am 2 - t) t-2M z 
Dropping corrections of order t/m 2 = 0(p 2 ), this expression simplifies to 
6(t-AM 2 ) f 2mVt-AM 2 2 



lmi(t) = IGnmVt \ arCtan t-2M 2 + ° {P } j ' (6) 
The problem addressed above shows up in this formula: The quantity 



_ 2mVt-AM 2 
X ~ t-2M 2 

represents a term of 0(l/p). The standard chiral expansion of InTy(t) thus 
corresponds to the series arctan x = |vr — l/x + |(l/x) 3 + ... This series, 
however, only converges for \x\ > 1. In the vicinity of t = 4M 2 , the condition 
is not met, so that the chiral expansion diverges. The problem arises because 
the quantity x takes small values there, while the low energy expansion treats 
it as a large term of 0(l/p). In the region |x| < 1, we may instead use the 
convergent series arctan x = x — | x 3 + . . . , but this amounts to an expansion 
in inverse powers of p. 

The rapid variation of the form factor near t = 4M 2 is related to the fact 
that the function arctan z exhibits branch points at z — ±i. The analytic 
continuation of j(t) to the second sheet therefore contains a branch point 
just below threshold: 

(t-AM 2 )(Am 2 -t) %r9 M 4 

= -1 -> t = AM 2 . 

(t - 2M 2 ) 2 m 2 

This implies that, in the threshold region, the form factor does not admit an 
expansion in powers of meson momenta and quark masses. As was shown 
in ref. the heavy baryon perturbation series to 0(p 3 ) coincides with the 
chiral expansion of the relativistic result and it was noted in ref. that 
this representation does not make sense near t = AM 2 . The corresponding 
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imaginary part amounts to the approximation arctan x — > ^tt, so that the 
singularity structure on the second sheet is discarded. Within HBXPT, an 
infinite series of internal line insertions must be summed up to properly 
describe the behaviour of the form factor near threshold. The relativistic 
formula (||), on the other hand, does apply in the entire low energy region, 
because it involves the full function j(t) rather than the first one or two 
terms in the chiral expansion thereof. 



4 Low energy expansion near threshold 

We conclude that in the threshold region, the low energy structure cannot 
be analyzed in terms of the standard chiral expansion. The first two terms 
of this expansion, i.e. the first two terms of the heavy baryon perturbation 
series, 

, . 6U-AM 2 ) (irM a(t-2M 2 )} 



Vi ^t{t-AM 2 )\ 



provide a decent representation only if t is not close to threshold. 

To resolve the structure in the threshold region, we need to consider an 
expansion that does not treat the quantity x as large. This can be done by 
replacing the variable t with the dimensionless parameter 

^ V t - AM 2 M 

* = u — ' « = — , 

a M m 
and expanding the amplitude at fixed £, so that the momentum transfer 

t = (4 + a 2 i 2 ) M 2 

stays close to threshold. 

For the imaginary part, the expansion at fixed £ takes the form 

T . . 6(t-4M 2 ) , . ^. 2N1 

Im7(t) = j 32^^ arcta ^ + ° (P) ^ 
At large values of £, this representation smoothly joins the one provided by 
the heavy baryon expansion, where t/M 2 is kept fixed. To amalgamate the 
two, we note that, in the threshold region, the quantity Im7i(t) reduces to 
the first two terms of the series arctan £ = | n — arctan l/£ = 1 7r — l/£ + . . . 
Hence the difference between Inry(i) and the first two terms of the heavy 
baryon perturbation series is approximately given by 

, . 6{t-AM 2 ) f olM aM 
Im7 2 (t) = — — < — — arctan ■ 



32nmM [ y/ 1 - AM 2 Vt — AM 2 
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Outside the threshold region, 111172 (t) is negligibly small - in the chiral count- 
ing of powers, it represents a term of order p 2 . The representation 

Im 7 (t) = Im 7l (t) + Im 72 (t) + 0(p) (9) 

holds irrespective of whether t/M 2 or £ is held fixed. Indeed, one may verify 
that on the entire interval 4M 2 < t < 20 M 2 , the formula differs from the 
expression in eq. (||) by less than 1%. 

Finally, we turn to the function 7 (t) itself. As mentioned above, the 
loop integral cannot be expressed in terms of elementary functions. We may, 
however, give an explicit representation that holds to first nonleading order of 
the low energy expansion, throughout the low energy region. The calculation 
is described in appendix A. It leads to a representation of the form 

7 (t)=7i(t)+72(*) + 0(p). (10) 

The first term corresponds to the result obtained in HBXPT ||. This piece 
explodes in the vicinity of the threshold, like its imaginary part. The explicit 
expression for the remainder reads 

1 \ olM f aM )) 

l2 (t) = <^ - In { 1 + I I . 11 

ny) 327rmM \ V4M 2 - 1 \ V4M 2 - 1 J J 1 ' 

In the language of HBXPT, this term is generated by multiple internal line 
insertions. It takes significant values only in the immediate vicinity of thresh- 
old, where it cures the deficiencies of 71 (t). Indeed, this term does account 
for the branch cut at t = AM 2 — M^/m 2 , which is missing in 71 (t). 



5 Infrared singularities in the self energy 

Having established the need to sum certain graphs of the heavy baryon chi- 
ral perturbation series to all orders, we now formulate a general method that 
leads to a representation where the relevant graphs are automatically ac- 
counted for. The method relies on dimensional regularization: We analyze 
the infrared singularities of the loop integrals for an arbitrary value of the 
dimension d. 

To explain the essence of the method, we first consider the simplest ex- 
ample: the self-energy graph shown in fig. |2|. We again focus on the corre- 
sponding scalar loop integral 

„ 1 r d d k 1 
H = — 



i J (2vr) d (M 2 - k 2 - it) (m 2 - (P - k) 2 - it) 
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Figure 2: Self energy 
5.1 Singular and regular parts 

The integral H = H(P 2 ,M 2 ,m 2 ) converges for d < 4. We need to analyze 
it for nucleon momenta close to the mass shell: P = mv + r, where v is a 
timelike unit vector and r is a small quantity of order p. It is convenient to 
work with the dimensionless variables 

P 2_ m 2_ M 2 M 
tt = , a= — , 12 

which represent quantities of order Q = 0(p°) and a = 0(p), respectively 
(recall that m is the larger one of the two masses, a ~ |). 

In the limit M — > 0, the integral develops an infrared singularity, gener- 
ated by small values of the variable of integration, k — 0{p). In that region, 
the first factor in the denominator is of 0(p 2 ), while the second is of order 
p. Accordingly, the chiral expansion of H contains contributions of O^^ -3 ). 
We may enhance these by considering small dimensions. For d < 3, the 
leading term in the chiral expansion of H exclusively stems from the region 
k = 0(p), which generates a singular contribution of order p d ~ 3 , as well as 
nonleading terms of order p d ~ 2 , p d_1 , . . . The remainder of the integration 
region does not contain infrared singularities and thus yields a contribution 
that can be expanded in an ordinary power series. For sufficiently large, 
negative values of d, the infrared region dominates the chiral expansion to 
any desired order. 

To work out the infrared singular piece, we use the standard Schwinger- 
Feynman-parametrization 



ab h {(1 



dz 



(13) 



z) a + z b } 2 
Performing the integration over k, we obtain 

H = KfdzCi- 2 , k= (47r)-3 m d - 4 r(2-f) , (14) 
Jo 

C = z 2 - 2 a Q z ( 1 - z) + a 2 (l - z) 2 - ie . 
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In this representation, the infrared singularity arises from small values of 
z: There, the factor C vanishes if a tends to zero. We may isolate the 
divergent part by scaling the variable of integration, z = au. The upper 
limit then becomes large. We extend the integration to oo and define the 
infrared singular part of the loop integral by 

/*oo poo 

I = k dzCi~ 2 = M " duD-2- 2 , (15) 
Jo Jo 

D = 1 - 2Qu + u 2 + 2au (Qu - 1) + a 2 u 2 - ie . 

For short, we also refer to / as the infrared part of H . The remainder is 
given by 

/oo 
dzCi~ 2 . (16) 

The decomposition^ 

H = I + R (17) 

neatly separates the infrared singular part from the regular part: For nonin- 
teger values of the dimension, the chiral expansion of / exclusively contains 
fractional powers of p, 

I = 0(p d ~ 3 ) + 0( P d ~ 2 ) + o( P d - 1 ) + . . . , 
while the corresponding expansion of R is an ordinary Taylor series, 

R = 0(p°) + 0(p 1 ) + 0(p 2 ) + ... 

5.2 Properties of the decomposition 

For an arbitrary value of d, the explicit expressions for the quantities H, I, 
R involve hypergeometric functions. In four dimensions, the corresponding 
integrals are elementary. We will give the explicit representations for d = 4 
when we discuss renormalization (see sections [| and|J). 

In terms of the dimensionless variables a and Q, the chiral expansion of 
the infrared part takes the form 

I = m d - 4 a d - 3 + at 1 + a 2 t 2 + ... ) . (18) 

3 In the terminology of Ellis and Tang Q , / represents the soft component of the ampli- 
tude, while R is the hard component. We do not use these terms to avoid confusion with 
the standard concepts, which concern the behaviour at short rather than long distances. 
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The coefficients i n , which only depend on f2, are obtained by expanding the 
integrand in eq. fll5l) in powers of a. The corresponding explicit expressions 
also involve hypergeometric functions. The expansion coefficients of the re- 
mainder, on the other hand, are simple polynomials of Q: 

m"r(2-j) f l + (rf-6)tf 1 

(te)l(rf-3) 1 (i-5) + "7 ' 1 ' 

The series contains poles at d = 3, 4, 5 . . . We now wish to show that, for 
d < 3, the chiral expansion of R converges throughout the low energy region. 

The integrand of the representation (|T6|) is analytic in a, except for the 
cuts associated with the zeros of C. These are located at 

a = -^—^ ± VQ 2 - l) . 
The expansion of the integrand thus converges in the disk 



a 



< | Q ± Vtt 2 - 1 1 . (20) 
The condition is obeyed for all real values of in the interval 

1 + a 2 „ 1 + a 2 

< Q < 



2a 2a 

This range corresponds to < P 2 < 2m 2 + 2M 2 and thus generously covers 
the entire low energy region. Moreover, for d < 3, the integral converges 



uniformly, so that R is analytic in a: The series (0) converges for all values 
of Q in the above interval. 

In the case of /, the representation (|lfj) is relevant. The zeros of the term 
D are located at 



a 



Q ± Vtt 2 - 1 . 



In the range — 1 < Q < 1, the square root is imaginary, so that the zeros 
occur at \a\ 2 = (1/u — Q) 2 + 1 — Q 2 . This expression has a minimum at 
u = On the entire interval of integration, the chiral expansion of the 

integrand thus converges if \a\ 2 < 1 — Q 2 , or, equivalently 



- VI -a 2 < n < VI - a 2 . (21) 

Again, the integral converges uniformly for d < 3. Hence the series ( |TB| ) 
converges if is in the range (^) . This interval is considerably more narrow 
than the one found above, but in the present context, the above result suffices: 
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It demonstrates that the two parts occurring in the decomposition of the self 
energy can unambiguously be characterized by their analytic properties at 
low energies. Since the functions /, R are analytic in the external momenta, 
their values are uniquely determined also outside the above region. Indeed, 
we will show in section || that the chiral expansion of / also converges in the 
entire low energy region. 



5.3 Representation in terms of modified propagators 

It is instructive to interpret the above decomposition in terms of the formula 
(Tl3|). The infrared part results if the integral is taken from to oo rather 
than from to 1, 

dz 1 



{ (1 — z)a + zb} 2 a(b — a) 



This shows that the decomposition (|T7| ) corresponds to the two terms in the 
algebraic identity 

~ h = ~7JT — V " ' ^ 

ab a[b — a) b[b — a) 

Note, however, that the ie-prescription drops out in the difference b — a, so 
that the integral over k of the individual terms on the right is ambiguous. To 
avoid the ambiguity, we may for instance equip the two masses with different 
imaginary parts. The proper expression for the infrared part reads 

1 f d d k 1 

(23) 



% J (2vr) d (M 2 - k 2 - ie) {m 2 - P 2 + 2Pk — M 2 — ie) ' 

It differs from H in that the term k 2 occurring in the nucleon propagator is 
replaced by M 2 . The regular part, on the other hand, is given by 

o = _l i m) 

% J (2vr) d (m 2 - (P - k) 2 - ie) (m 2 - P 2 + 2Pk — M 2 — ie) K ' 

The expression involves two heavy particle propagators, one where the term 
with k 2 is retained, one where this term is replaced by M 2 . 

The function H is symmetric with respect to an interchange of M and 
m. The term / collects those contributions of the expansion in powers of 
a = Mj m that involve fractional powers of M, while R collects the fractional 
powers of m. The operation M ^ m, k ^ P — k interchanges a and b, so 
that the two terms on the r.h.s. of eq. (0) are mapped into one another. 
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This suggests that an interchange of the two masses takes / into R and vice 
versa. The argument does not go through, however: The difference in the ie- 
prescriptions shows that the operation does not take the integral in eq. ( |2"3"| ) 
into the one in eq. (0). The sum / + R, of course, is symmetric under an 
interchange of the two masses. 

In the heavy baryon chiral perturbation series, the scalar self energy graph 
of fig. |2| is replaced by an infinite string of one loop graphs, involving an 

/ \ => / \ + / \ + / \ + ... 
-* — > ^- — X — )( )( — 

Figure 3: Internal line insertions. The double line denotes the heavy baryon 
propagator [2mv ■ (k — r) — it}' 1 , the cross an insertion of (k — r) 2 . 

arbitrary number of internal line insertions (fig. [|). The corresponding scalar 
loop integrals are obtained from the relativistic version with P = mv + r, 
treating both r and k as small quantities of order p and expanding the nucleon 
propagator in powers of p: 



2mv ■ {k — r) — {k — r) 2 — ie 



1 (k-r) 2 

+ tt. -T-, s t^t + 



2mv ■ (k — r) — ie {2mv ■ (k — r) — ie} 2 

The integral over the leading term converges for d < 3, yielding a contribution 
of order p d ~ 3 that depends on the vector r only through the projection v ■ r. 
The integral over the second term of the series converges for d < 2 and yields 
a term of order p d ~ 2 , etc. The individual contributions are not Lorentz 
invariant, but the series may be reordered in such a manner that only the 
Lorentz invariant combination 2m v ■ r + r 2 = P 2 — m 2 enters: The heavy 
baryon series then reproduces the expansion of I in eq. (|T8|), term by term. 

To demonstrate that this is so, we recall that the infrared part dominates 
the chiral expansion to any desired order if d is taken sufficiently negative: For 
d < 3 — n, the region k = 0(p) yields all of the terms in the chiral expansion 
of the integral H, up to and including p d ~ 3 + n . In that region of integration, 
however, it is legitimate to interchange the integration with the expansion. 
This is precisely what is done in the heavy baryon approach. Hence that 
approach does yield the expansion of J, to any finite order: The infrared 
part of the relativistic loop integral represents the sum of the corresponding 
integrals occurring in the heavy baryon series. The difference between the 
two formulations of the effective theory resides in the regular part: In the 
heavy baryon approach, this part is absent. The representation shows 
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that R corresponds to a loop formed with two nucleon lines - evaluating 
integrals of this type in the manner described above, the rules of HBXPT 
]T5| indeed yield R = 0, order by order. 



5.4 Dispersive representation 

The function H(P 2 ,M 2 ,m 2 ) obeys a dispersion relation, which in four di- 
mensions requires one subtraction. Setting s = P 2 and suppressing the other 
variables, the relation takes the form 

H{s) = H[so) + _± | (3 ,_ 3o)(3 ,_ 3 _ 16) ^W) , (25) 
Im H(s) = e ( s ~ s +) , s± = (m± M) 2 . 

16 7T S 

The function p(s) stands for the familiar two-body phase space factor 



p(s) = y/(s-s+)(s-s_) = 2 MmVtt 2 - 1 . 

The logarithmic divergence of the loop integral manifests itself in the sub- 
traction constant H(sq), which contains a pole at d = 4. If we subtract at 
threshold, sq = s+, the subtraction constant is given 



r(2-f) m ^-3 + M d - 3 

~ (47r)l(d-3) (m + M) ' 

The structure of this expression is typical: It contains a term with a fractional 
power of M and one with a fractional power of m, representing the infrared 
singular and regular parts, respectively, 

r(2-|) m^- 3 



(47r)i(d-3) (m + M) ' 

r(2-|) 



(47r)l(d-3) (m + M) ' 

The dispersion relation ( p5| ) is a variant of the formula ( ]I4]) : In the limit 
d — >• 4, the corresponding representation for H(s) — H(s+) is proportional to 
dz\n(C/C + ), where C + is the value of C at s = s + . Since C is linear in s, 



C7 = ^—^-( s '- s -ie), 

m 2 M 2 ds' \dz\ 

s = + 



4 While the representation for the subtraction term holds in any dimension, the one for 
the discontinuity is valid only in the limit d — > 4. 
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the argument of the logarithm reduces to (s' — s — ie)/{s' — s + ). With an 
integration by parts this indeed leads to eq. (|25D . Note that, on the interval 
< z < 1, the function s' = s'(z) has a minimum at 



M 



m + M 



s' ■ = s 

"mm 



The map z — > s' thus covers the interval (m + M) 2 < s' < oo twice. 

In the Feynman parameter representation, the only difference between 
H, I and R is that the integrations extend over different intervals. The 
one relevant for R is 1 < z < oo, which is mapped onto — oo < s' < 0. 
Accordingly, R has a cut along the negative real axis, but is analytic in the 
right half plane. The discontinuity across the cut is only half as big as in the 
case of H, because the interval is now covered only once: 

ImR(s) P{S) - - 



32 7TS 



The expression only holds for s / 0: In addition to the cut, the function 
R(s) contains a pole at s — 0. In view of the singular behaviour of the 
discontinuity, the dispersion relation cannot be written in the form (p5|). 
Instead, we may establish a twice subtracted dispersion relation for s R(s) 
(compare eq. ( p?T| ) in section |8|): 

s + m 2 -M 2 , (s-s + )(s-s-) y° ds' 



s 32 7r 2 s J _ oa p(s') (s' — s — it) 

m + M , , 

In contrast to the subtraction term in the dispersion relation for the full 
integral H(s), the one occurring here does not represent a constant, but 
contains a pole at s = 0. That point, however, is far outside the region 
covered by chiral perturbation theory. Neither the subtraction term nor 
the dispersion integral contain singularities in the low energy region. The 
dispersive representation neatly demonstrates that the chiral expansion of R 
is an ordinary Taylor series. 

The function / is the difference between H and R and hence has a cut 
on the left as well as one on the right: 

1mI{a) = ^L Us-s + )-h{-a) 
lbrcs I 2 

Note that the infrared part possesses the same discontinuity across the right 
hand cut as the full integral, even far away from threshold. 
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6 Generalization 



We now generalize the above analysis to arbitrary one loop graphs. All of 
these can be reduced to integrals of the form 



i J (2rr) d a 1 ...a m b 1 ...b n ' 
The denominator stems from the meson and nucleon propagators: 

a, = M 2 -{k- qi f - ie , 

b, j = m 2 - (Pj - k) 2 - ie . 

Part of the numerator is generated by the derivative couplings characteristic 
of chiral perturbation theory. The remainder arises from the term ft that 
occurs in the numerator of the nucleon propagator, Pj —I/t+m. The external 
meson momenta q? represent quantities of order p. The nucleon momenta 
P? are close to the mass shell, Pj = m? + 0(p). 

On account of Lorentz invariance, the above integral may be decomposed 
in a basis formed with tensor polynomials of the external momenta P/*, q^. 
Inverting this representation, the coefficients of the decomposition may be 
expressed in terms of scalar integrals, where the momentum factors are re- 
placed by their projections onto the external momenta or by factors of k 2 . 
These, however, represent linear combinations of the terms a«, bk occurring 
in the denominator:^ 



If the graph in question involves several different external momenta, the 
above procedure may be rather clumsy, but it shows that all one loop integrals 
arising in the 7riV-system may be expressed in terms of the scalar functions 



6.1 Singular and regular parts 

The self energy corresponds to H = Hn, the triangle graph to 7 = H 2 i. 
Concerning the counting of powers of momentum, loops that do not involve 

5 Without loss of generality, we may put one of the external momenta to zero, say q±. 
6 A more efficient method is described in ref. Jig . As shown there, the tensor integrals 
may be generated by applying suitable differential operators to the scalar ones. 




M 2 — di , P x -k = \{bx-a x + Pi -m 2 + M 2 ) , 



mn i J (2n) d a\ . . . a m b\ . . . b, 
so that it suffices to study the properties of these.Q 




(26) 
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the propagation of a heavy particle are trivial: The integral then represents a 
homogeneous function of order that a regular part does not occur, 

ImO = H m0 , R m0 = . 

In the opposite extreme, where the graph exclusively involves nucleon 
propagators, we may shift the variable of integration with k! = k — P\. 
The denominator then depends on the external momenta only through the 
differences Pj — Pi, which represent small terms of 0(p). In Euclidean space, 
the integrand is thus approximately given by (m 2 + k' 2 )~ n , so that the region 
k' = 0(m) dominates. For d < 2n—r, the first r terms of the chiral expansion 
may be worked out by performing the expansion under the integral - the 
coefficients are polynomials in the external momenta. Hence the integral 
does not contain any infrared singularities, 



Ion — , Ron — H{ 



On 



If the loop contains meson as well as nucleon propagators, the integral 
involves both an infrared singular and a regular piece. The infrared singu- 
larities arise from the region k = 0{p). There, each of the pion propagators 
yields a factor of p~ 2 , while each of the nucleon propagators yields a factor 
of p~ l , so that the infrared region generates a contribution of order p d ~ 2m - n . 

We may combine all of the mesonic propagators by means of the formula 

d \ (m_1) f\ f\ X 

dxi . . . dx m _i 



a x ...a m \ dM 2 J J " ' J A 

The numerator is given by (m > 2 and X = 1 if m = 2) 

X = x 2 ■ (x 3 ) 2 ■ ■ ■ (x m ^) m ' 2 . 

The denominator is obtained recursively, with 

Ap+i = x p A p + (1 - x p ) a p+1 (p = 1, . . . , m - 1) 

Ai = ai , A = A m . 

The result for A is quadratic in k, 

A = A - (k - q) 2 - it . 

The constant term A is of order p 2 , while q represents a linear combination 
of external momenta and is of order p. Likewise, if there are several nucleon 
propagators, we may combine these with 

1 / d f 1 Y 
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In this case, the denominator is of the form 

B = B - (P - kf - ie , 
with P 2 = m 2 + 0(p), B = m 2 + 0(p). The loop integral then becomes 

H m =(~T^(~T 1) f 1 d^yXY 1 1 



dM 2 J V dm 2 J J tf i J (2n) d AB 



where we have abbreviated the integral over the m + n — 2 Feynman param- 
eters by L dxdy. The integral over k is the one studied in section 



1 [AJLJ_ = H{P ^M 2 ,m^ 
J (2n) d AB K ' 



i J (2TT) d AB ' \ P^P-q, M 2 ^A,m' 2 ^B ' 

so that the analysis given there may be taken over. The decomposition 
H = I + R leads to an analogous splitting for the general scalar one loop 
integral: 

In particular, the representation for I mn reads 

I mn = j^dxdyXY(- M ) l((P-q) 2 ,A,B) . (27) 

In view of the formula (P^D, the representation of the infrared part in terms 
of Feynman parameters coincides with the one obtained for the full integral 
H mn , except that the integration over one of these parameters - the one 
needed to combine the meson propagators with the nucleon propagators - 
runs from to oo rather than from to 1. Performing the derivatives with 
respect to the masses and integrating over k, we obtain 



T(m + n — ~ 



a\ pi poo 

^ / dxdyXY dz{\-z) m ~ x z n ~ x C%- 
Jo Jo 



(An) 

C = (1 - z) A + z B - z (1 - z) {P - q) 2 - it . (28) 

The analogous representation for R mn is obtained by replacing the interval 
of integration for z by 1 < z < oo and changing the overall sign. 
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6.2 Uniqueness of the decomposition 



Formally, we may expand this representation in the same manner as the self 
energy: Set Pj = mv + Tj, rescale the variable of integration with z = au 
and expand the integrand in powers of q^, rj and M. The same result is 
obtained, term by term, if we treat the loop momentum k in the integral (|26|) 
as a quantity of 0(p) and perform the chiral expansion under the integral. 
Applying the procedure also to R mn , we obtain two series of the form 

Iran — TTl Ct ( i mn + OL% mn + U i mn + . . . ^ , (29) 

i? = m d ~ 2m ~ 2n ( r° A-ar 1 + a 2 r 2 + I 

1 *>mn nL I 1 ran ~ " ' mn ' u ' mn ~ ■ ■ ■ J ■ 

While the dimensionless coefficients i v mn are nontrivial functions of the vari- 
ables qi/M, . . . , q m /M and ri/M, . . . , r n /M, those occurring in the expan- 
sion of Rmn are polynomials. 

As discussed in detail in section |3], the chiral expansion in general makes 
sense only in part of the low energy region. In the representation (EST), the 



problem arises from the fact that it is not always legitimate to interchange the 
integration over the Feynman parameters with this expansion - the integral 
I mn as such describes the low energy behaviour perfectly well. 

We again emphasize that the decomposition of the integral into an in- 
frared part and a remainder is uniquely characterized by the analytic prop- 
erties of the two pieces, also if the expansion (|29"D only holds in part of the 
low energy region. The proof closely follows the one given in section 5.2 and 
we only indicate the modifications needed to adapt it to the present more 
general situation. As shown there, the chiral expansion of / and R does 
converge in part of the low energy region. In the present context, the range 
corresponds to 

- 2 {A (B - A)} 1 < (P - q) 2 - A - B < 2 {A (B - A)} 1 . (30) 

The quantities A, B, q and P depend on the parameters x\, ... 

... , y n -i used to combine the meson and nucleon propagators, respec- 



tively - the condition (p0|) must be met for all values of these parameters 
in the range < Xi,yj < 1. If this is the case, the chiral expansion of the 
integrand on the r.h.s. of eq. (|27f) converges. 

It is not difficult to see that the low energy region does contain a domain 
where the above condition is met. The vector q is contained in the polyhedron 
spanned by the corners qi and analogously for P. We may, for instance, take 
all of the qi to be much smaller than M and put the in the immediate 
vicinity of a vector P that sits on the nucleon mass shell, such that A ~ M 2 , 
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B ~ m 2 , (P — g) 2 ~ m 2 . The condition is then obviously met in the entire 
region of integration. Moreover, the integral over the Feynman parameters 
converges uniformly there. This completes the proof. 

7 Comparison with HBXPT 

In the framework of the relativistic effective theory, the evaluation of an 
amplitude to one loop thus yields three categories of contributions, arising 
from (a) tree graphs, (b) infrared singular part and (c) regular part of the 
one loop integrals. At a given order of the chiral expansion, the one particle 
irreducible components of the regular parts are polynomials in the exter- 
nal momenta. In coordinate space, these contributions thus represent local 
terms: They are equivalent to the tree graph contributions generated by a 
suitable Lagrangian, AC. So, if we replace the effective Lagrangian by 

we may drop the regular parts of the loop integrals. The resulting represen- 
tation is identical to the one obtained in the heavy baryon approach, except 
that the infrared parts of the one loop graphs are included to all orders - 
the problems afflicting the chiral expansion of the infrared singularities are 
avoided. We add a few remarks concerning the above relation between the 
term C e ff relevant for the original form of the effective theory and the effective 
Lagrangian occurring in our framework, C' e ^. 

7.1 Lorentz invariance 

First, we note that both of these schemes are characterized by a Lorentz 
invariant effective Lagrangian: By construction, the term AC is Lorentz in- 
variant. In explicit formulations of HBXPT, the invariance of the effective 
Lagrangian is by no means manifest. One of the reasons is that the equa- 
tions of motion are used to eliminate two of the four components of the Dirac 
spinor that describes the nucleon in the relativistic formulation of the the- 
ory. In fact, it is perfectly legitimate to use the equations of motion: The 
resulting modification of the effective Lagrangian is equivalent to a change 
of variables. The operation, however, destroys manifest Lorentz invariance 
- in terms of the new variables, the transformation law of the field takes a 
rather complicated form. 

The point here is that all of that can be avoided. Instead of explicitly 
performing the chiral expansion of the Lagrangian and evaluating the per- 
turbation series with the corresponding nonrelativistic propagators, we may 
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simply replace the integrands of the various loop integrals by their chiral 
series, i.e. perform the nonrelativistic expansion before doing the integral [§. 
As discussed above, this procedure amounts to replacing the relativistic loop 
integrals by their infrared parts. The result for the various amplitudes of in- 
terest is the same as the one obtained within the standard approach, except 
that our method accounts for the mass insertions to all orders. 

7.2 Chiral symmetry 

Both the relativistic and the heavy baryon formulations of the effective theory 
are based on an effective Lagrangian that is manifestly invariant under chiral 
transformations. In the above construction, this property of the term C' e ^ 
is not evident. Although the equivalence with the standard heavy baryon 
approach ensures chiral symmetry, it is instructive to see how this property 
arises within the present framework. 

For this purpose, we first formulate chiral symmetry in terms of objects 
that are amenable to an evaluation in the framework of the effective the- 
ory: Consider all Green functions of the type (N(P r , s')\T 0\ . . . O r \N(P, s)), 
where the operators Oi represent vector, axial, scalar or pseudoscalar quark 
currents. Chiral symmetry implies that, in the limit where the quark masses 
are put equal to zero, these matrix elements are interrelated through a set of 
Ward identities. We now analyze the implications of these identities for the 
regular parts of the one loop graphs responsible for the term AC. 

At one loop, the Green functions are represented by a sum of three contri- 
butions belonging to the three categories a, b, c specified above. Since £ e // 
is invariant under chiral symmetry, the tree graph contributions (a) obey 
the Ward identities. Furthermore, dimensional regularization preserves the 
symmetries of the Lagrangian. Hence the contributions from the one loop 
graphs (b + c) also obey these identities, for any value of the regularization 
parameter d. Now, the chiral expansion of the infrared singular (regular) part 
only contains fractional (integer) powers of the chiral expansion parameter 
p. Hence the Ward identities can be satisfied by the sum of the two pieces 
only if they are obeyed separately by the two parts: The vertices collected in 
AC obey the same set of linear constraints as the vertices contained in £ e ff- 
This explains why the term AC is invariant under chiral transformations^. 

7 For the case of the purely mesonic vertices, it is explicitly demonstrated in ref. |l7| 
that the Ward identities indeed imply a symmetric effective Lagrangian, but we did not 
perform the corresponding analysis for the present, more general case. 
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7.3 Infrared part as an alternative regularization 



Since C e ff contains all terms permitted by Lorentz invariance and chiral sym- 
metry, the modification £ e ff — > £ e ff + AC is equivalent to a change of the 
effective coupling constants: The bare coupling constants to be used in the 
original form of the relativistic effective theory differ from those occurring 
in our scheme. In this respect, the two methods of calculation appear like 
two different regularizations of the theory - for the physical amplitudes to 
be independent thereof, the values of the bare couplings must be tuned to 
the regularization used. 

Note, however, that the two prescriptions for the evaluation of the loop 
integrals in general lead to different results even if these are convergent. View- 
ing H^ n -^ r and as two different regularizations of the same integral, 
we are leaving the standard class of admissible regularizations (dimensional, 
Pauli-Villars, momentum space cutoff, . . . ). 

For a Lagrangian that contains all of the Lorentz invariant vertices that 
can be formed with the field and its derivatives, it would be natural to take 
the loop integrals as being defined only modulo an arbitrary Lorentz invari- 
ant polynomial. In our context, that class is too large. Since the effective 
Lagrangian is chirally invariant, the same coupling constant determines the 
strength of an entire string of vertices. A change in one of the couplings 
generates a specific polynomial contribution in several different amplitudes. 
Conversely, if we are removing a polynomial from one of the loop integrals, 
we need to subtract a corresponding term in some of the other loop integrals, 
too - otherwise, the procedure would yield amplitudes that do not obey all 
of the Ward identities. It is essential here that dimensional regularization 
preserves the symmetry and that this procedure allows us to unambiguously 
identify the infrared part for all of the integrals. At any finite order of the chi- 
ral expansion, the difference between the dimensional regularization H^ n -^ r 
and the infrared part Z^"^ is a polynomial and the polynomials occurring 
in different loop integrals are correlated in such a manner that the Ward 
identities are obeyed: If we consistently replace all of the integrals by their 
infrared parts, the content of the theory remains the same. Hence, it is le- 
gitimate to think of the infrared part as an alternative regularization of the 
loop integrals and to indicate this with the symbol J T d d k : 

Mr _ 1 r d d k k^ ... k^ 
mn = ~i JjJZkY a l ...a m b l ...b n ' 

The standard regularizations yield the smoothest possible high energy 
behaviour: The expansion in inverse powers of p starts with p d+r ~ 2m ~ 2n . The 
high energy behaviour is crucial for renormalizability but in the context of 
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effective low energy theories, it is irrelevant, because the high energy domain 
is anyway outside the reach of the framework. The infrared regularization 
Imri'^ r instead yields maximally smooth behaviour at low energies: All of 
the regular contributions of order p°, p 1 , ... are absorbed in the effective 
coupling constants, so that the low energy expansion starts with an infrared 
singular piece of order p d + r - 2m - n _ 

In principle, the analysis given in the preceding section may be extended 
to arbitrary graphs. The leading infrared singularity originates in the region 
where all of the loop momenta are small. Disregarding momentum factors 
that may arise from the vertices or from the propagators, the leading infrared 
singularity occurring in the low energy expansion of a graph with I loops, 
m mesonic and n baryonic propagators is of order p £d - 2m - n - the counting 
of powers is the same as in heavy baryon chiral perturbation theory. At the 
present stage of our understanding, the extension beyond the one loop ap- 
proximation is a rather academic issue, however. In the case of 7rN scattering 
for instance, significant progress could be achieved by extending the known 
results to order p A . With the method outlined above, this should require 
rather little effort: It suffices to (i) replace the dimensional regularization 
used in ref. JI| by the infrared regularization and (ii) add the contributions 
from the one loop graphs generated by £ 2 . 



8 Renormalization 

We now consider the behaviour of the loop integrals in the limit d — > 4 and 
again start with the self energy. 



8.1 Self energy 

The integral over the Feynman parameter z in eq. (|T5| ) only converges for 
d < 3. The continuation to d — > 4 may be performed as follows. The factor 
C is of the form 

C = Cq + C\(z — Zq) 2 , Zq — 



l + 2«n + a 2 ' 
C ^TTT^?- ^( 1 + 2aW). 

Replace C^~ 2 by CoC%~ 3 + C\(z — zo) 2 C^~ 3 . The second term is propor- 
tional to the derivative of C*2~ 2 . An integration by parts leads to 

■ d -2 _ K J „.d-4 _ , I j A\n / 



k\ dzCi- z = - — - { a d - A z + (d - 4) C / dzC*- 6 } . (31) 



" ° L. ^0 
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Since the remaining integral converges for d < 5, the right hand side can now 
be continued analytically to d — 4. The factor k contains a pole there, 



K=£ iy^- 4 =- 2 *-!^ +o(d - 4) - 
x -s{^-^" ite+r(i)+i 

We have expressed the singularity in terms of the standard pole term A, 
which contains a running scale fi. The scale relevant here is the mass of the 
nucleon: A represents the value of A at the scale /j, — m. The renormalized 
amplitude, which we denote by /, is obtained by removing the pole 

I = I - \ u , 

P 2 - m 2 + M 2 



P 2 

For the regular part and for the full integral, the renormalizations read 

R = R + (2 + v) A , H = H + 2X. 
By construction, we have 

H =I+R. 



The counter term for H is momentum independent, but the quantity v is 
not. As v does not contain infrared singularities, its expansion in powers of 
M 2 and P 2 — m 2 is an ordinary Taylor series. 



8.2 Renormalization scale 

It is important here that, in the relativistic formulation of the effective the- 
ory, loops involving nucleon propagators contain an intrinsic scale, even in 
the chiral limit: the nucleon mass. This is in marked contrast to the mesonic 
sector and to the standard heavy baryon approach, where dimensionally reg- 
ularized loop integrals are scale invariant in the chiral limit, so that the 
removal of the divergences necessarily involves a free parameter, the running 
scale. As is well known, this does not give rise to ambiguities, because the 
renormalization of the loop integrals only requires polynomial counter terms: 
It suffices to tune the coupling constants of those terms in the effective La- 
grangian that enter at the order of the chiral expansion considered - the 
result for quantities of physical interest then becomes scale independent. 
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In the present context, the situation is different. The chiral expansion 
of the infrared part I contains arbitrarily high orders - in the language of 
standard HBXPT, we are summing up an infinite number of graphs. Their 
renormalization requires counter term polynomials of arbitrarily high order: 
Heavy baryon graphs of 0(p n ) call for counter terms of 0(p n ). Indeed, the 
chiral expansion of the counter term A v contains polynomials of arbitrarily 
high order. If we were to work with an arbitrary running scale, we would 
need to include infinitely many terms in the effective Lagrangian and tune 
their scale dependence properly - only then the amplitudes would become 
scale independent. 

Neither can this be done in practice, nor is it necessary. All of the loop 
integrals that require counter terms with a nonpolynomial momentum depen- 
dence contain an intrinsic scale and we may identify the renormalization scale 
with this one, i.e. set fi = m. A running scale is needed only for loops formed 
exclusively with mesonic propagators - these do not contain an intrinsic one. 

8.3 Other one loop graphs 

The generalization to other one loop graphs is straightforward. Concerning 
the full scalar integrals, only the self energy H requires renormalization: 
For m + n > 2 the functions H mn represent convergent integrals in four 
dimensions. Nevertheless, the infrared parts of H12 and H21 do contain a 
pole at d — 4, because the integral over z in the representation (|28|) only 
converges for m + n > 3. We define the renormalized infrared parts by 

For m + n > 2, the full integral converges, so that 

Rmn Rmn A v mn if m + n > 2 . 

As shown in section |5], H mn , I mn and R mn may be represented as integrals 
over a derivative of H, I and R, respectively. The counter terms for I mn or 
Rmn are thus given by a derivative of v with respect to the masses. In accord 
with the statements made above, v is linear in M 2 and m 2 , so that 

v mn = if m > 2 or n > 2 . 

In the case of I12 (triangle graph with one meson and two nucleon propaga- 
tors), the formula ( ^7|) involves a single Feynman parameter: 

I12 = ~ J dy^((P-q)\A,B) , 

Ql = q = Q, P = yP 1 + (l-y)P 2 , 

A = M 2 , B = m 2 - y(l - y)(P 1 - P 2 ) 2 . 
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The counter term is thus given by 



^12 



V 



dy 1 arctanh 77 



{y Pl + (l-y )P 2 } 2 ~ (P 1 -P 2 ) T) 



As was to be expected on general grounds, the counter term does not con- 
tain any infrared singularities. Since r\ represents a term of 0(p), the chiral 
expansion starts with 

i ' 12 = -(i^) +0(!,2, ■ 

The renormalization of the function J 2 i (triangle with two meson and 
one nucleon propagators) may be worked out in the same manner. The 
corresponding counter term is given by 

This completes the list of renormalizations for the scalar loop integrals. 



9 Convergence of the chiral expansion, 
explicit representations 

In section |5], we have shown that, in the case of the self energy, the chiral 
expansion of the infrared part converges if the variable Q is in the range 
(21). We now wish to show that the expansion actually converges throughout 



the low energy region. This is most easily done on the basis of an explicit 
representation. 

9.1 Explicit representation of the self energy 

At d = 4, the integral remaining on the r.h.s. of eq. (|3T| ) is elementary. In 
terms of the variables Q and a of eq. (|T2|), the result for the renormalized 
infrared part reads (— 1 <fl < 1) 

- 1 aVl-W ( Q + a \ 

arccos I -^^^=^^^= I (66) 



1 a{Q + a) (21na-l) 



a 



16tt 2 l + 2afi + a 2 
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In accord with the power counting of HBXPT, the chiral expansion of I starts 
at order p. The coefficients are nontrivial functions of Q: 



16tt 2 



{2 Vl - ^ 2 arccos(-fi) +Q(2\na- 1)J + 0( 



a 2 ) 



For the regular part and for the full integral, the explicit expressions may be 
written in the form ( — 1 <0< 1): 



- 1 a^T^Q 2 . ( a^T^Q 2 

R = — arcsm 



+ 



n 2 l + 2an + a 2 VVl + 2aVL + a 2 

1 1 + aQ 



16vr 2 1 + 2 afi + a 



2 - 



1 aVl-O 2 

il = ^ — arccosl— iz) 

8tt 2 1 + 2 a + a 2 v ; 

1 a(fz + a) , 1 
In a + 



3vr 2 l + 2«n + a 2 167T 2 



In agreement with eq. (|T9|), the chiral expansion of the regular part starts at 
0(p°) and only contains polynomials, 

R= -^-{1-aQ + a 2 } + 0(a 3 ) . 

l07T z 



9.2 Chiral expansion of the self energy 

Consider now the expansion of the function H in powers of a at fixed Q. The 
radius of convergence is determined by the zeros of the denominator, which 
occur at a = —Q ± \/fi 2 — 1. Hence the chiral expansion of H converges in 
the disk ( p0|) - the convergence region is the same as the one relevant for 
R. In view of I = H — R, the statement also holds for the expansion of the 
infrared part. This proves the claim made above. 

It is essential here that we consider the chiral expansion at fixed Q. If we 
instead set P = mv + r, expand in powers of r and M and collect terms of 
the same order, a phenomenon similar to the one encountered in the chiral 
expansion of 7(i) occurs. The ordering of the double series amounts to setting 
r = Mr and expanding in powers of M at fixed f. The parametrization 
implies 

Q = v . f + — (f 2 - 1) . 
2m 

In contrast to the expansion considered above, where Q stays put, we are 
now expanding this variable around the value v ■ f. 
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The loop integral has a branch point at threshold, P 2 = (m + M) 2 . In 
the variable Q, this singularity occurs at O = 1. In the limit M — > 0, the 
branch point is mapped into the plane v ■ f — 1. Accordingly, the radius 
of convergence becomes small if f happens to be close to this plane: The 
coefficients of the expansion blow up if v ■ f tends to 1. In other words, the 
series only converges in part of the low energy region. This illustrates the 
fact that the convergence of the nonrelativistic expansion is a rather delicate 
matter, sensitive to the details of the infrared structure. 

9.3 Other one loop integrals 

As shown in section ||2, there is a range of external momenta where the chiral 
expansion converges, for all one loop graphs. The example of the triangle 
graph shows, however, that the low energy region in general contains holes, 
where the chiral expansion breaks down. Even for the self energy, we need 
to order the series suitably for the expansion to converge throughout the low 
energy region. In that case, this can be done by expanding at fixed fl We 
do not know of a corresponding set of variables for the general loop integral. 

Our method does not rely on the chiral expansion of the loop integrals - 
on the contrary, the problems afflicting that expansion motivated the present 
work. We have reformulated HBXPT in such a manner that the relevant 
infrared singularities are summed up. 

For loop integrals with more than two vertices, the explicit representation 
becomes complicated. As is well known, all of the one loop graphs can be 
expressed in terms of dilogarithms. A much simpler representation may, 
however, be given if we resort to the approximation discussed in section |], 
which amounts to summing up only the leading infrared singularities. There, 
we studied the low energy properties of the function ^(t) - in the above 
terminology, this function coincides with H21, except that the two external 
nucleon momenta are put on the mass shell (P 2 = m 2 , q 2 = t, Pq = \ t). 
In that case, we considered two different expansions, one at fixed t/M 2 , the 
other at fixed (t — AM^/M 4 and then joined the two. For a corresponding 
representation of the infrared part, we refer to the appendix. An alternative 
procedure might be to look for a uniformizing variable replacing t. The 
breakdown of the chiral expansion is generated by the fact that the form 
factor contains a branch point both on the first and on the second sheet and 
that the two move together if the mass of the nucleon is sent to infinity. 

As shown in appendix [FJ, the sum of the leading infrared singularities of 
all one loop integrals that are relevant for the scalar form factor and for the 
elastic nN scattering amplitude can be represented in terms of elementary 
functions, but we cannot offer a general method that would lead to such a 
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representation. In a given case, the issue may be settled by trial and error, for 
instance by numerically comparing the infrared part in the kinematic region 
of interest with the first one or two terms in the chiral expansion thereof. If 
the expansion fails, we may search for an improved approximation - this is 
what we did in the case of the triangle graph. 



10 Chiral expansion of the nucleon mass 

As a first illustration of our method, we now evaluate the physical mass of 
the nucleon to order p 4 . The two-point-function of the field ip(x) may be 
represented in the form 

i/d d xe lPx (0|T^(x)^(0)|0) ' 



m 



The leading contribution to E = S(P) is of order p 2 - it stems from the 

,(2) 

'AT 



term 4 ciM 2 ipijj contained in C^- The one loop graphs shown in fig. are 



generated by C$ , and start contributing at order p 3 . Finally, there is a 
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Figure 4: One loop graphs contributing to the self energy of the nucleon. 
The cross denotes a vertex from C$ . 



tree graph contribution from C$ (see eq. (0)): 

S = -4 Cl M 2 + S a + S b + S c + eiM 4 + 0(p 5 ) . 

The explicit expressions for the loop contributions are obtained with the 
standard rules of relativistic perturbation theory: 

S a = ^k( m + f){M 2 I+(m-f)flW-A N } , 

3M 2 A 7T f_ P 2 
F 2 

E c = -4ciM 2 



r2 5S a 



dm 
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The mass insertion in graph c arises from the shift m — ► m 2 = m — Ac\M 2 in 

(2) 

the nucleon propagator, generated by the term from C N mentioned above. 
We could have replaced the mass in the free part of the Lagrangian by m 2 , 
so that the graph f|c would then be absent. 

The only difference to the standard evaluation of the graphs is that the 
loop integrals are regularized in a different manner. In particular, the full 
scalar self energy integral (fig. §a) is replaced by the infrared part thereof, 
/ = I(P 2 ,M 2 ,m 2 ) - we have discussed the properties of this function in 
detail in the preceding sections. The term A^v denotes the scalar nucleon 
propagator at the origin, 



A 



1 r d d k 1 



N — — 



i J j (2ir) d m 2 — k 2 — it 

In infrared regularization, that term vanishes, because it does not contain 
any infrared singularities (see section ^ : 

A N = . 

For purely mesonic loops such as the one occurring in fig. ^b, there is no dif- 
ference between dimensional and infrared regularization: The integral J 7 d d k 
coincides with the ordinary ci-dimensional integral. The graph is proportional 
to the pion propagator at the origin, which we denote by A^, 

d d k 1 n% ,2 / T , 1 , M\ 



A n = - - — -t — = = 2M Z X H In — . 

i J j (2n) d M 2 -k 2 -ie \ 16n 2 m J 

Finally, the integral I^ 1 ' may be expressed in terms of / (see appendix [B]) : 
/ (1) = ^ {(P 2 -m 2 + M 2 )/ + A.-A^} . 

The same relation also holds in dimensional regularization. 

The physical mass of the nucleon, which we denote by m^r, is determined 
by the position of the pole in the two-point-function. In view of the factor 
m — f 1 , the term 1^ does not contribute. Evaluating the quantity / with 
the explicit expression in eq. (0), we obtain 

m N = m _4 Cl M 2 - yA \ +fc 1 M 4 ln — + fc 2 M 4 + 0(M 5 ) , 

32ttF 2 m 

3 

kl = ~ on 2P2 (9a - § Cim + c 2 m + 4 c 3 m) , (34) 
3 

k2 = El ~ i no 2F2 ( 2 9 A ~ C 2 m ) " 
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The formula agrees with the result of refs. fll|, |12||. The bare constants M, F, 
m, gA, ci, c 2 , c 3 , c 4 remain finite when the regularization is removed, but e\ 
contains a pole at d — 4. The quantity ei is the corresponding renormalized 
coupling constant at scale fi = m : 

3A 

e~\ = e x - (V4 - 8 cim + c 2 m + 4 c 3 m) . (35) 

11 Wave function renormalization 

The wave function renormalization constant is the residue of the pole in the 
two-point-function and is determined by a derivative of the self energy with 
respect to the momentum, 



=m N 



With the above expression for E, which is valid to order p 4 , we can extract 
the residue to accuracy p 3 . The result, 

1 LM 1 *M\l +0(Mr 



2F 2 { 16tt 2 V m 3 2m, , 

is in agreement with those obtained within the heavy baryon formalism. For 
a detailed discussion of the latter, see ref. |~2] . 



Note that the multiplicative renormalization i\f en = Z~^i\) does not ren- 
der the two-point-function finite at d — 4. The reason is the following. We 
may collect all of the correlation functions associated with ip and ip by adding 
a term of the form fj ip + ip rj to the effective Lagrangian, where r](x) is an 
anticommuting external field. Such a term, however, breaks chiral symmetry: 
Under a chiral rotation, the field ip transforms with a matrix that involves 
the pion field. We cannot subject rj to the same rotation, because this field 
stays put when the meson variables are integrated out. So, off the mass 
shell, the correlation functions of the effective field are regularization depen- 
dent objects. In the context of the effective field theory, these are without 
significance - the field ip{x) merely represents a variable of integration. 

Instead we could consider the field \1/ = uipR + vft tpL, which does trans- 
form with a factor that is independent of the meson field: ^ — > Vr 

—>■ Vl^l- Accordingly, the counter terms needed to renormalize the 
quantum fluctuations generated by the term fj ^ + ^ rj are chirally invariant. 
In contrast to the correlation functions of ip, those of \l/ can unambiguously 
be worked out. Indeed, these variables are relevant for the low energy anal- 
ysis of QCD operators that are formed with three quark fields and carry 
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the quantum numbers of the nucleon: The leading term in the effective field 
theory representation of such an operator is a multiple of 

The main point is that the correlation functions of ip do not represent 
physical quantities. The graphs for say the scalar form factor or the nucleon 
mass also involve nucleons propagating off the mass shell, but the result is 
uniquely determined by the Lagrangian: The low energy representation of 
these quantities in terms of the renormalized coupling constants is indepen- 
dent of the regularization used. The one for the correlation functions of the 
effective field is unambiguous only on the mass shell. On-shell matrix el- 
ements such as form factors or scattering amplitudes may be obtained by 
extracting the residues of the relevant poles in suitable correlation functions 
- in that connection, only the on-shell properties matter. 

12 Scalar form factor to order p A 

In section [|, we considered the low energy representation of the scalar form 
factor to 0(p 3 ). We now extend that representation to the next order of 
the expansion (for an analogous calculation of the scalar form factors within 
SU(3), based on HBXPT, see ref. [@). For this purpose, we treat the quark 
masses m u , as external fields and calculate the response of the transition 
amplitude (N(P', s') out|iV(P, s) in) to a local change in these fields. Within 
the effective theory, the transition amplitude may be worked out by treating 
the term x in the effective Lagrangian as a space-time dependent quantity and 
evaluating the two-point-function (0 out|T-?/>(:r) ip(y) \ in) in the presence of 
this external field: The transition amplitude is determined by the residue 
of the double pole occurring in the Fourier transform of this quantity at 
P 2 = m 2 N , P' 2 = m 2 N . To extract the residue, we amputate the external 
nucleon legs, evaluate the remainder at P 2 = P' 2 = m 2 N and multiply the 
result with the wave function renormalization constant Z. 

12.1 Evaluation of the graphs 

The tree graph contributions from and read: 

a tree = Z (-4ciM 2 + 2 dM 4 + e 2 M 2 t) . 

The one loop graphs are shown in fig. [5[ All of these, except the one in 
fig. ||a represent contributions of 0(p 4 ). When performing the perturbative 
calculation, it is convenient to include the term 4 ciM 2 ip ip in the free part of 
the Lagrangian, replacing m by m 2 = m — Ac\M 2 . The nucleons occurring in 
the various graphs then propagate with m 2 . Note that the difference m 2 — P 2 
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Figure 5: The loop graphs contributing to the scalar form factor at 0(p A ). 



represents a quantity of order p. In the propagator, the shift m — > m 2 thus 
generates a first order correction. The distinction between m, and itln 
only matters in the triangle diagram [5|a. In fact, graph [5]b represents the 
change occurring in this diagram if m is replaced by m-i - this graph is absent 
if the mass insertion from c\ is included in the free part of the Lagrangian. 
In view of rriN — m2 = 0(p 3 ) we may also replace the factors of arising 
from the matrix elements u'fu and u'f'u by m 2 . The contributions from 
the one loop diagrams then take the form: 



3M 2 g 2 A 



Or 



2F 2 

6ciM 4 T , , 3c 2 M 2 



{4m s lg\t)-mJ(t)\ 



F 2 v ' 8F 2 m 2 
3c 3 M 2 



IF 2 



{(t-2M 2 ) J(t) + 2A W } 



6ciM 2 



3Cl9 p A 2 M {4M 2 m 2 I l2 (t) - Am 2 I^(m 2 ) + A.} 



The notation is specified in appendix |B[ With the relations given there, 
the various invariants may be expressed in terms of the basic functions 
J(t), I(s), hiit) and In(t). The individual terms entering the combina- 
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tion 4m 3 /g ) (t) -mJ(t) are of order p°, but the leading terms cancel: The 
low energy expansion of a a only starts at 0(p 3 ), in accord with the counting 
of powers in graph |5]a. Hence we may ignore the difference between m2 and 
m in the above expression. 

The loop graphs contain divergences proportional to M 4 and to M 2 t, 
respectively. The renormalization ( |35l) of the coupling constant e\ removes 
the first one. The second one requires the following renormalization of e 2 : 

e 2 = e 2 + or i {3 g\ + c 2 m + 6 c 3 m} A . 
2b A m 

12.2 Result 

We write the result for a(t) = a tree + a a + + cr c + + a e in the form 

a{t) = a + a{t) , 

where a = cr(0) represents the value of the form factor at the origin and 
is referred to as the cr-term. According to the Feynman-Hellmann theorem, 
this term represents the derivative of the nucleon mass with respect to the 
quark masses m u , rrid, or, equivalently, 

Indeed, the sum of the contributions from the various graphs agrees with the 
derivative of the formula ([34]) for the nucleon mass. 
For the remainder, the calculation yields: 

a{t) = 3 ^ 2 2m {(t - 2M 2 ) I 21 (t) + 2M 2 / 21 (0)} (36) 
+{k 3 M A + k 4 M 2 t} J(t) + k 5 MH + 0(p 5 ) , 
h = --^2 i 6c i - c 2 - 3c 3 } , 

k 4 = - — - {3 g\ + 2 c 2 m + 12 c 3 m} , 

1 1 
fc s = ~^fc 4 lna+ {9g 2 A + 2c 2 m + 36 c 3 m} + e 2 . 

The only difference to the result of an analogous calculation in HBXPT is 
that the representation for the function I 2 i(t) given in appendix || also covers 
the vicinity of t = 4M 2 , where the heavy baryon representation fails. The 
value of (x(t) at the Cheng-Dashen point t = 2M 2 was given earlier, in ref. (19 
- our expression confirms this result. 
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12.3 Unitarity 

Unitarity offers an instructive test of the momentum dependence. Within the 
effective theory, the contributions from intermediate states containing more 
than two pions only show up at three loop order. Hence the representation 



obtained for the form factor must obey the elastic unitarity relation |20 



lma(t) = 9(t - 4 Ml) 3 f a*(t) f + (t) + 0(p 7 ) , (37) 

where cr n (t) is the form factor associated with the cr-term of the pion, 

a w (t) = (7T (p')\m{uu + dd)\7c {p)) 

and fl(t) is the t-channel I = J = partial wave amplitude of irN scattering. 

To calculate the left hand side of the unitarity condition, we recall that 
the function l2i(t) represents the renormalized infrared part of the triangle 
integral j(t) introduced in section [3[ At any finite order of the low energy 
expansion, the difference between the infrared part and the full integral is a 
polynomial, so that Im/ 2 i(^) = hajit). The explicit expression was given in 
eq. (H). For J(t), the imaginary part reads 



lmJ(t) = 9(t-4M-) \^ ■ 

The quantities on the right hand side of the unitarity condition are needed 
only at tree level, where a n (t) = M 2 . The corresponding approximation 
for the ttN scattering amplitude is given in appendix |C[ The comparison 
shows that the representation obtained for the form factor indeed obeys the 
unitarity condition, up to contributions that are beyond the accuracy of a 
one loop calculation. 

12.4 Value at the Cheng-Dashen point 

The low energy theorem that underlies determinations of the cx-term from 
ttN data relates the scattering amplitude to the scalar form factor at the 
Cheng-Dashen point, where t = 2M 2 . It is therefore of interest to evaluate 
the difference 

A CT = a(2M 2 J - a(0) 
with the above representation of the form factor. The result is of the form 

A a = Ai M 3 + A 2 M 4 ln — + A 3 M 4 + 0(M 5 ) , 

m 

A 3 ^ 
1 64ttF 2 ' 
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The terms of order M 4 involve the coupling constants of £ 



(2). 
N ■ 



A _ 3 9A C 2 3C 3 /oo\ 

2 " 167r 2 F 2 m 16vr 2 F 2 8tt 2 F 2 ' 1 ' 

3(2 + tt)^ 3(4-7r) Cl (14-37r)c 2 3c 3 

3 1287r 2 F 2 m 16tt 2 F 2 192tt 2 F 2 16tt 2 F 2 62 ' 

To express these in terms of observable quantities, we compare the tree graphs 
for the ttN scattering amplitude generated by + with the subthresh- 
old expansion of Holder and collaborators |2l| (see appendix 0). Since that 



comparison allows us to determine the effective coupling constants only up 
to corrections of order M, we denote the resulting estimates by c$ : 



4M2 «o + 2M 2 d+ } , 4 U) = ^ d+ (39) 



JO) _ ^ / ,+ , on,f2 >+l 



c (o) _ _ F 2 ,+ (o) _ J_ r F 2 _ 11 

c 3 — r "01 J c 4 — 2 m I 00 2 J ' 

The above expressions for A 2 and A3 then take the form 
a _ 3g A d 10 ^doi 



l67T 2 F 2 m 32tt 2 8tt 2 ' 
3(2 + 7r)^ 2 3(4-vr)4 (14 - 3vr) d+ 3(2 - tt) d+ 
3 128vr 2 F 2 m 64vr 2 M 2 384vr 2 32tt 2 2 ' 

The difference between F, F.,, and M, is beyond the accuracy of the 
representation (^). We replace F by F n = 92.4 MeV and use the mass of 
the charged pion. For the same reason, we may identify the bare coupling 
constant g A with the experimental value g A = 1.267 ± 0.0035 |2"2| . Using the 
values of the subthreshold coefficients quoted in appendix O, we then obtain 
Ax M 3 = 7.6 MeV, A 2 M 4 In a = 7.8 MeV and (A 3 - 2 e 2 ) M 4 = -1.4 MeV. 
These terms add up to 

A,,- = 14.0 MeV + 2 M 4 e 2 , 



to be compared with the result of the dispersive calculation of ref. [20] : 



A a = 15.2 ± 0.4 MeV . (40) 

The comparison shows that the contribution from the coupling constant e 2 
is small, as it should be. 

The above calculation resolves an old puzzle: The leading term in the 
chiral expansion of A CT - the one of order M 3 - accounts for only half of the 
result. The terms of order M 4 are numerically of the same size, because they 
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Figure 6: The t-channel I = J = partial wave The normalization 

factor n is chosen such that A CT is the area under the curve (dotted: Born 
approximation, dash-dotted: tree graphs of , full: dispersion theory, 

all quantities in units of M n ). 



are enhanced by a chiral logarithm. The underlying physics can be sorted out 
by noting that the elastic unitarity condition (|3"7|) leads to the representation 



A a = / dtn(t)o*(t)f°(t) , (41) 

J4M2 



3M 2 y/t-AM* 

n(t) — 



Tit? {t-2M 2 ) (Am 2 N -t) 

The one loop calculation discussed above replaces the pion a-term by the 
first term in the chiral series, cr^(t) — > M 2 . The contribution proportional 
to g\ represents the value of the dispersion integral that results if the par- 
tial wave amplitude is replaced by the Born approximation fg + (t) = 0(p), 
given in eq. ( |C.2| ). The remainder is dominated by the coupling constant C3, 
which generates a polynomial contribution^ to the partial wave amplitude: 
f°(t) C3 = (2M 2 - t) C3 m 2 / AttF 2 = 0(p 2 ). 

The tree approximation for /?(£) is shown in fig. |^, together with the 

Born term (tree graph from alone) . To indicate the weight in the disper- 
sion integral (f4~ID, we plot the corresponding contributions to the quantity 

8 Since the corresponding contribution to the scattering amplitude only depends on t, 
the t-channel partial wave projection is trivial and is obtained by multiplying the relevant 
term in D + (v,t) with the factor (4m 2 - i)/16vr = m 2 /4vr + 0{p 2 ). 
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n(t) f+(t). We also depict the result of the dispersive evaluation described in 
ref. pOj , which includes the higher orders of the chiral series^ and is therefore 
complex - the curve shown represents the quantity n(t) Qualitatively, 
the picture is quite similar to the one found for the imaginary part of the 
scalar form factor of the pion, Imovft): The higher order effects tend to 



amplify the leading order terms also in that case |23| , [24 . 

The figure demonstrates that, at low energies, the tree approximation to 
the scattering amplitude provides a rather decent representation. The Born 
term alone, however, only dominates in the immediate vicinity of t = 4M^, 
where it exhibits the peculiar structure discussed in section [|. For \fi > 
380 MeV, the contribution generated by the coupling constants ci, c 2 and C3 
is more important than the one proportional to g\. 

Although the straightforward expansion of A a in powers of the quark 
masses is well-defined and can unambiguously be worked out, the first term 
of that series does not yield a decent approximation. The term arises from 
the infrared singularity generated by the Born approximation. On the one 
hand, the contributions from the second term of the chiral expansion are 
suppressed by one power of M, on the other, they are enhanced by a chiral 
logarithm - numerically, they are equally important. 

12.5 The (T-term 

Finally, we turn to the value of the form factor at t = 0. The first four 
terms in the chiral expansion of the nucleon mass were given in eq. fl34"|). The 
corresponding expansion for a = er(0) is of the form 

a = aiM 2 + a 2 M 3 + 03M 4 In — + a 4 M 4 + 0(M 5 ) . (42) 

m 

The coefficients follow from the Feynman-Hellmann theorem: 

a-i = -4 Cl , (43) 
9 Si 



a 2 



64vrF 2 



3 

°3 = -1 p 2F2 — (g\-Bc 1 m + c 2 m + Ac 3 m) , 
3 

°"4 = - rA 9F9 — (3 9a~ 8cim + 4c 3 m) + 2ei 



The tree approximation (39) for the effective coupling constants Ci, c 2 , C3 
suffices to work out the numerical values of a 2 , cr 3 , as well as the corre- 
sponding contribution to 04. With the values of the subthreshold coefficients 

9 The value quoted in eq. (^0|) also accounts for the higher order terms in the form factor 
ov(i), which we are ignoring here, because they start showing up in A CT only at 0(p 5 ). 
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quoted in appendix |C], we obtain a 2 M 3 = -23MeV, a 3 M 4 \na = -9 MeV, 
(04 — 2ei)M 4 = 0.9 MeV. To evaluate 0\ in an analogous manner, we 
need a more accurate determination of ci, but we may instead take the 



experimental value a = 45 MeV from ref . |20| . This leads to the estimate 
crxM 2 + 2eiM 4 = 76MeV, subject to an uncertainty of about 8 MeV. 

These numbers show that the expansion of a in powers of the quark 
masses contains a large contribution from the infrared singularity generated 
by the Born term, 02 M 3 . At the next order of the expansion, there is a 
logarithmic infrared singularity, cr 3 M 4 In a. The comparison with the leading 
term o\M 2 shows that this effect is of the size typical for chiral logarithms. 
The main difference to the situation encountered in the mesonic sector is 
that the expansion contains odd as well as even powers of p. If the series 
is truncated at order p 3 , we must expect a less accurate representation than 
the one obtained in the mesonic sector at one loop. 

The determination of the coupling constant c\ to the accuracy needed 
here requires a calculation of the ttN scattering amplitude to order p 4 and is 
beyond the scope of the present paper. We can, however, study the effects 
generated by the leading infrared singularities in the values of the effective 
coupling constants. For this purpose, we make use of the results reported in 
ref. |1|, ||, 0, where the scattering amplitude is evaluated to order p 3 . The 
formulae for the coefficients of the subthreshold expansion given in ref. |7], P5fl 
imply that the corrections of 0(p 3 ) generate the following shifts in the values 
of the effective coupling constants: 

cf ) = c f ) + 6 Ci (44) 

1536vrF 2 ' 64vrF 2 

5c 3 = - — — — - M , <5c 4 = A A - M . 
6 768ttF 2 ' 16vrF 2 

As a check, we have applied the method discussed in section || to the loop 
integrals that occur in the representation of the scattering amplitude given 
in ref. [[[]]. In the vicinity of the point v — t — 0, the chiral expansion 
converges for all of these (see appendix 0). The result confirms the formulae 
for the coefficients of the subthreshold expansion quoted above and thus 
also corroborates the expressions fl44|) for the first order shifts in the effective 
coupling constants. The one for c\ allows us to establish the relation between 
the er-term and the subthreshold coefficients to accuracy p 3 : 



a = Fl (4 + 2MX) + ll^l + 0{Mf, 
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to be compared with the prediction at tree level, a = F% (g?q + 2M^g?q~ 1 ), 
where the third term is missing. Numerically, the shift amounts to only 
2.1 MeV: In the relation between the observables a, g?qo and d^, the effects 
generated by the leading infrared singularities are an order of magnitude 
smaller than those seen in the chiral expansion of a. This demonstrates that 
the large infrared singular contribution that occurs in the chiral expansion 
of the u-term at 0(M 3 ) also appears in the scattering amplitude - the 
relevant combination of subthreshold coefficients picks up nearly the same 
contribution. Note that the magnitude of the a-term as such is not at issue 
here. For the values of the coefficients d^ , d^ given in ref. jn[], the tree level 
and one loop results for a are 50 and 52 MeV, respectively. We could have 
used a somewhat different input, for instance one for which the outcome 
for the cr-term is 45 MeV - the difference between the tree and one loop 
approximations would be the same. 



13 Discussion 

13.1 Reordering of the perturbation series 

The fact that the bulk of the infrared singular contributions occurring in the 
chiral expansion of the nucleon mass and of the cr-term merely amounts to a 
shift of the bare parameters was noted long ago EfJ. In the preceding section, 
we listed the change in the effective coupling constants ci, C2, C3, C4 generated 
by the infrared singularities at 0(p 3 ) and pointed out that, in the relations 
between the observable quantities considered, these singularities nearly can- 
cel. We expect a similar, albeit less dramatic reduction also to occur in other 
quantities of physical interest. Some of the fluctuations seen in the published 
results for the threshold parameters of ttN scattering, for instance, merely 
reflect the fact that the expansion of the various observables in powers of the 
quark masses contain large contributions from infrared singularities. 

A more significant comparison of the results obtained at tree level and 
at one loop results if one considers the predictions of the effective theory: 
Instead of comparing the contributions arising in the chiral expansion of a 
given observable at various orders of the expansion, one may compare the 
relations between observable quantities which follow at tree level with those 
obtained at one loop. The effective coupling constants entering the one loop 
result are different from those relevant at tree level. The tree level values 
are the best choice for the representation of the scattering amplitude to 
0(p 2 ), those in eq. ([H]) are relevant for the improved representation of this 
amplitude that accounts for the terms of 0(p 3 ). 
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The elimination of the coupling constants in favour of physical quantities 
amounts to a reordering of the perturbation series. In more explicit terms, 
the reordering we are advocating here is the following (we restrict ourselves to 
a discussion to one loop). The perturbation series relies on a decomposition 
of the effective Lagrangian into a leading term and a perturbation, 



The standard ordering results if Co is identified with the first two terms of 
the derivative expansion. We argue that the decomposition 



with Set = Ci — q , is a better choice, because it reduces the magnitude of 
the perturbation C\. Like for the standard bookkeeping, an evaluation of 
the various observables to 0{p 2 ) requires the calculation of the tree graphs 
belonging to Co; for a representation to 0(p 4 ), we need to add the tree graphs 
of Ci as well as the one loop graphs of Cq. 

Two modifications occur when the perturbation series is extended from 
0{p 2 ) to 0(p 4 ): new graphs and the change 5ci in the values of the coupling 
constants (the explicit expressions for the latter in eq. ([H]) account for the 
shifts of order 5ci = O(M), but do not include those of order M 2 ). As 
discussed above, the two effects partly cancel. In the above bookkeeping, 
the change 5ci is booked among the corrections, on the same footing as the 
contributions from the loop graphs, which both affect the observables and 
are responsible for the shift in the coupling constants. The cancellations 
between the two types of contributions thus tend to reduce the magnitude 
of the perturbations generated by C%. 

At first sight, the claim that the values of the coupling constants depend 
on the order at which the perturbation series is considered, may appear 
to contradict the fact that these constants represent perfectly well-defined 
quantities that determine the chiral expansion coefficients of the various ob- 
servables. As is well known, the results found for these coefficients on the 
basis of a calculation to some given order of chiral perturbation theory repre- 
sent low energy theorems that remain strictly valid if the expansion is carried 
to higher orders. We do not put this into question, but merely emphasize 
that the result obtained for the values of the coupling constants does depend 
on the order to which the perturbation series is worked out. 

To illustrate the need for a distinction between the coupling constants as 
such and their values at a given order of the perturbation series, we again 
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consider the quantity A a and focus on the chiral logarithm contained therein, 
A 2 M 4 lna. The term arises from the infrared singularity occurring at the 
lower end of the dispersion integral (pEIp. According to eq. Q58D, the coefficient 
A 2 contains a piece that is proportional to the coupling constants c 2 and c 3 . 
It arises from the contribution to the partial wave /+(£) that is generated 

by the tree graphs of C N . Hence the coupling constants relevant for an 
evaluation of A 2 are those pertaining to the tree level representation of the 
scattering amplitude. If we were to calculate A 2 with the improved values 
of the coupling constants that follow from the relations (f0|), we would in 
effect be using a tree level representation for /?(£) with the wrong coupling 
constants. The point here is that the same one loop graphs that give rise to 
a change in the values of c 2 , C3 also modify the partial wave amplitude. Since 
a substantial fraction of the infrared singularities occurring at 0(p 3 ) affects 
the term A a and the scattering amplitude in the same way, it does not make 
sense to account for one of the changes without accounting for the other. 

In a certain sense, the odd and even powers of the chiral series lead 
a life of their own. The quantity A CT illustrates the fact that the leading 
terms in both of the subseries need to be investigated in order to arrive at a 
significant result. To our knowledge, the available one loop results for the ttN 
scattering amplitude account for the expansion only to 0(p 3 ). It will be very 
interesting to compare the predictions that follow from the representation 
of the amplitude to 0(p 4 ) with those obtained from the same experimental 
input at tree level. 

The problem also occurs in the mesonic sector. The infrared singularities 
are weaker there, because the expansion only involves even powers of p. At 
the precision reached with a two loop calculation, however, the need for a 
distinction between the coupling constants as such and the values found at 
a given order of the perturbation series manifests itself quite clearly. As 
discussed in refs. [^7|, |28|, for instance, inconsistencies may arise if the two 
loop representation for the 7T7t scattering amplitude is evaluated with the 
values of the coupling constants £1, £ 2 obtained on the basis of a one loop 
calculation from K ei decay. 

13.2 The role of the A(1232) 

The occurrence of a chiral logarithm explains why an evaluation of A CT to 
order M 3 does not yield a decent approximation for this quantity. The coef- 
ficient of the logarithm involves the combination c 2 + 6 C3 of effective coupling 
constants, which is dominated by C3. It is understood why the value of this 
coupling constant obtained from low energy ttN phenomenology is large: This 
constant receives an important contribution from the singularities generated 
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by the A(1232). In fact, it was noted in ref. |29] that a calculation of the 
scalar form factor that explicitly includes the A degrees of freedom yields a 
result for A CT that is consistent with the one obtained from dispersion theory. 
The role of this state for the low energy structure in the baryonic sector 



is discussed in detail in the literature || |2T], |29|, 0, |3T| . For a recent review, 
in particular also of the small scale expansion that allows one to analyze 
the extension of the effective theory in a controlled manner, we refer to |32| . 



The extended theory is compared with the framework used in the preceding 
sections in appendix [D], where we also estimate the contributions to the 
effective coupling constants that are generated by the A. 

For the quantities analyzed in the present paper, it is not essential whether 
the A is incorporated as a dynamical variable in the effective Lagrangian or 
whether the effects generated by this state are accounted for only indirectly, 
through the values of the effective coupling constants: In the domain studied 
here, the graphs describing the exchange of a A are adequately described by 
those terms of the chiral expansion that occur up to and including 0(p 4 ). 

In the Mandelstam plane, the point around which the chiral expansion 
is performed corresponds to Sq = m 2 + M 2 , t = 0. The expansion of the 
resonance denominators is controlled by the ratio x = (s — So) 2 /( m A — s o) 2 - 
For small values of s — Sq and t, in particular near the Adler zero and for 
quantities like cr(0) or cr(2M 2 ), the expansion rapidly converges. At the 
threshold, where x = 4M 2 m 2 /(m\ — m 2 — M 2 ) 2 ~ 0.18, the expansion of 
the resonance denominators is still under good control. For higher energies, 
however, the singularities generated by the A must explicitly be accounted 
for to arrive at a decent representation of the scattering amplitude. 

In the mesonic sector, the p plays an analogous role. Remarkably, as far as 
the expansion of the resonance denominators is concerned, the convergence 
radius is roughly the same: m\ — m 2 — M 2 ~ m 2 — 2M 2 (the left hand 
side is even a little larger). The two Mandelstam triangles, however, are 
of very different size: The one relevant for ittt scattering is smaller by the 
factor M/m, so that the threshold is much closer to the Adler zero. At the 
threshold, the parameter that controls the expansion of the term l/(m 2 — s) 
is very small indeed: x = 4M 4 /(m 2 - 2M 2 ) 2 ~ 0.005. In the mesonic sector, 
an effective theory that does not explicitly account for this singularity yields 
meaningful results even well above threshold. 

13.3 Comparison with the static model 

It is instructive to compare the framework discussed above with the earliest 
version of an effective field theory for the baryons, the static model. This 
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model is characterized by the Hamiltonian [ |33|j : 

H = j d\^[-k 2 + V7T 2 + M ff 2 7T 2 } + Jd 3 X^p{x)^{(T-V) (TT.T)V. 

The nucleon is kept fixed at the origin, ?/>(a;) — > ip(t). The corresponding 
four-dimensional space of states is spanned by vectors that differ in the spin 
and isospin quantum numbers - the operator ip < . . . tp generates transitions 
between these. The result for the self energy of the nucleon reads 



E = jd\d\ Vp(f) • Vp{y) - 



-M v \x—y\ 



x-y\ 

The function p(x) is normalized to J d 3 x p(x) = 1 - it describes the structure 
of a nucleon that is stripped of its meson cloud. 

In the relativistic formulation of the effective theory, the nucleon is rep- 
resented as a point particle, p(x) — > 5(x), and the leading contributions to 
the self energy of the nucleon arise from one loop graphs. To compare the 
expressions for the self energy, we note that in the static model, the expan- 
sion of the exponential yields the series 1 — M n r + ^M 2 r 2 — |M^r 3 + . . . 
The second term vanishes upon integration, so that the series takes the form 
E = Ei + E 2 M% + E 3 M 3 + ... In our formulation of BXPT, only the infrared 
part of the loop integrals is retained, the remainder being absorbed in the 
effective coupling constants. In that language, the term E\ is included in 
the bare mass m, while E 2 represents a contribution to the effective coupling 
constant c\. The coefficient E 3 is independent of the shape of p(x) and the 
value E 3 = —3g\/32irF% indeed reproduces the coefficient occurring in the 
chiral expansion (|34]) of the nucleon mass. At the next order of the expansion, 
however, the static model fails: The self energy E does not contain a chiral 
logarithm. In fact, in the local limit p(x) — > 5(x), the expansion terminates 
at 0{Ml). The deficiency arises because the model only accounts for the 
leading term C$ in the derivative expansion of the effective Lagrangian - 

the formula ( [43|) shows that, in the static limit (m — > 00), the nucleon mass 

(2) 

does not contain a chiral logarithm if the coupling constants Ci,c 2 ,c 3 6 C N 
are turned off. 

Quite a few other phenomena (ttN scattering, magnetic moments and 
electromagnetic form factors, pion photoproduction, Compton scattering, 
nuclear forces, for instance) were analyzed in detail within the static model, 
which provides a simple intuitive picture for the basic low energy features 
(for an excellent overview, we refer to the book of Henley and Thirring [ftSfl). 
A comparison of the results obtained with the modern version of the effective 
theory with those found in this model would be most instructive. 
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13.4 Momentum space cutoff 



The formulation of the effective theory in terms of pointlike nucleons has 
the advantage of being model independent. It systematically accounts for all 
contributions arising to a given order in the chiral expansion of the various 
observables. The price to pay is that the loop integrals cannot be inter- 
preted directly in physical terms. The machinery does incorporate the finite 
extension of the nucleon, but only indirectly, through the effective coupling 
constants. A more intuitive picture of the meson cloud that does explicitly 
account for the finite size of the nucleon was proposed in ref. |2(J. In that 
framework, the pointlike vertices are replaced by form factors, which in ef- 
fect cut the virtual meson momenta off. The proposal is based on the static 
model, where the nN interaction is also equipped with a form factor, given by 
the Fourier transform of the function p(x). The motivation for introducing 
such a cutoff was that the matrix elements of the perturbations encountered 
in the extension of BXPT from SU(2) to SU(3) are large. It does not make 
much sense to truncate the straightforward expansion of the baryon masses 
in powers of m Ul md and m s at order p 3 . The net effect of a cutoff is quali- 
tatively similar to the one resulting from the reordering discussed above: It 
reduces the magnitude of those terms that are treated as perturbations. 
The proposal is taken up in the recent literature In particular, 



the quantity A CT is studied within that framework in ref. |pq| . In the lan- 
guage of the dispersion relation (fll]), the calculation reported there essen- 
tially amounts to replacing the partial wave amplitude /+(£) by the Born 
term (dotted line in fig. |5p and cutting the dispersion integral off at mod- 
erate values of t. Such a calculation evidently yields a much smaller value 
for A CT than the one that follows, either from dispersion theory or from the 
evaluation of the chiral perturbation series to order p 4 . The analysis of sec- 
tion 0.4 shows that the contributions of order p 3 and p 4 are of a different 
origin. The terms of order p 4 are not properly accounted for by cutting off 
the virtual momenta in the one loop graphs of order p 3 . 



More generally, the problem with the approach proposed in ref. [26] is 
that it is model dependent. Moreover, introducing a cutoff in general ruins 
the Ward identities of chiral symmetry. The model independent method 
we are proposing does preserve the Ward identities as well as the infrared 
structure. The reordering of the perturbation series also leads to a more 
rapidly convergent expansion. It yet remains to be seen, however, whether 
this method will provide a coherent understanding of the mass spectrum of 
the baryon octet. 
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14 Conclusion 



1. We have shown that BXPT can be formulated in such a manner that 
both Lorentz invariance and chiral power counting are preserved at every 
stage of the calculation. The method relies on the fact that - for noninteger 
values of the dimension - the infrared singular parts of the loop graphs can 
unambiguously be separated from the remainder: They involve fractional 
powers of the chiral expansion parameters, while the remainder admits an 
ordinary Taylor series. The two parts are chirally invariant by themselves in 
the sense that they separately obey the Ward identities of chiral symmetry. 
This allows us to retain only the infrared singular parts of the loop integrals, 
absorbing the remainder in the coupling constants of the effective Lagrangian. 

2. The calculations required by our method are nearly identical to those 
relevant within the relativistic formulation of the effective theory given in 
ref. 0. The Lagrangian and the Feynman graphs are the same. The only 
difference is that the loop integrals are replaced by the corresponding infrared 
parts. At one loop order, this is achieved by extending one of the Feynman 
parameter integrations from the interval < z < 1 to the interval < z < oo. 

3. We have shown that, in the sector with baryon number 1, the chiral 
expansion breaks down in certain regions of phase space. The nonrelativis- 
tic expansion that underlies HBXPT inherits this problem and also leads to 
technical difficulties of its own, related to the fact that the nonrelativistic ex- 
pansion of infrared singularities is a rather subtle matter. A coherent heavy 
baryon representation only results if the insertions required by relativistic 
kinematics are summed up to all orders. The problem arises from the inter- 
change of the loop integration with the nonrelativistic expansion, which is 
not always legitimate. The method we are proposing avoids these difficulties 
ab initio, because it does not rely on a nonrelativistic expansion of the loop 
integrals. 

4. One of the advantages of chiral perturbation theory is that, at the 
algebraic level, the results are unambiguous. This implies that the chiral 
expansion of the representations obtained with the method proposed here 
must agree with the results of HBXPT, order by order - even in those cases 
where this expansion does not converge. 

5. As our method avoids the nonrelativistic expansion of the vertices, the 
number of graphs to be evaluated is smaller than in HBXPT. The price to 
pay is that the simplifications offered by nonrelativistic kinematics cannot be 
made use of. As an example, we mention processes involving external pho- 
tons, such as pion photoproduction. Bernard, Kaiser, Kambor and Meissner 
have shown that, in the heavy baryon approach, many of the one loop 
graphs occurring in that context can be dropped: In the nonrelativistic limit 
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and in the Coulomb gauge, vertices involving the coupling of photons to 
a nucleon line are suppressed. Since our method leads to the same power 
counting rules as the nonrelativistic expansion, we expect these simplifica- 
tions to also apply within our framework, but a relativistic formulation of 
the corresponding selection rules yet remains to be given. 

6. In the baryonic sector, the expansion involves odd as well as even 
powers of momenta. Most of the available calculations only account for the 
terms arising up to and including 0(p 3 ). The corresponding representation 
for the observables of physical interest is inherently less accurate than the one 
obtained at one loop in the mesonic sector, which holds up to and including 
0(p 4 ). A full one loop calculation of the ttN scattering amplitude yet remains 
to be carried out. 

7. We have illustrated our method with an evaluation of the cr-term and of 
the corresponding form factor to order p 4 . The form factor can be expressed 
in terms of elementary functions throughout the low energy region. In the 
domain where the standard chiral expansion in powers of momenta and quark 
masses is convergent, the result agrees with the one obtained on the basis of 
HBXPT. The representation constructed in the present paper, however, also 
holds in the vicinity of t = 4 M% , where the heavy baryon chiral perturbation 
series diverges. 

8. The example of the form factor shows that the expansion of some of 
the observables in powers of the quark masses contains a large contribution 
from the infrared singularities generated by the Born term, proportional to 
g\M 3 . We have shown, however, that the bulk of this contribution is common 
to the cr-term and to the nN scattering amplitude, so that it drops out 
when considering the predictions of the theory, i.e. the relations between the 
observables. 

9. Although the straightforward expansion of the various observables in 
powers of momenta and quark masses is perfectly meaningful, it does not 
yield a suitable ordering of the perturbation series. We have proposed a 
model independent reordering of this series, which can also be performed 
in HBXPT. At least in the two cases considered in the present paper (cr- 
term and momentum dependence of the scalar form factor), the operation 
strongly reduces the magnitude of the corrections. To see whether or not this 
is a peculiarity of the observables considered here, it would be of interest to 
perform the operation with the known representation of the tcN scattering 
amplitude to 0(p 3 ) [g § g g. 

10. We have discussed the effects generated by the A(1232) in some 
detail. This state plays a role similar to the one of the p in the mesonic 
sector: In the Mandelstam plane, the singularities associated with the A 
and the p occur at about the same distance from the point around which 
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the chiral expansion is performed. As far as the behaviour of the scattering 
amplitude in the vicinity of the Adler zero is concerned, or for quantities 
such as <t(0) or <r(2M 2 ), the singularities generated by the A are adequately 
described by their contributions to the effective coupling constants. 

11. The physical region of ttN scattering is further away from the Adler 
zero than the physical region of irn scattering, by one power of m/M. Al- 
though, at the threshold, the singularities generated by the A may still be 
replaced by their contributions to the effective coupling constants, we must 
expect the perturbation series for the threshold parameters to converge less 
rapidly for nN scattering than for tttt scattering. 

12. In the experimentally accessible region, the A does play a prominent 
role. The work reported in ref. [^] indicates that, even in the vicinity of the 
resonance, a decent description of the observed behaviour of the scattering 
amplitude may be obtained by supplementing the tree graphs of a simple 
resonance model with unitarity corrections. It yet remains to be seen whether 
a systematic analysis of the extended effective theory that includes the A 
among the dynamical variables and incorporates chiral symmetry ab initio 
will allow us to establish direct contact with the wealth of experimental data. 

13. The problems discussed in the present paper are not peculiar to the 
nucleons, but occur whenever the effective theory contains degrees of freedom 
that remain massive in the chiral limit. An application of our method to 7rK 
scattering within SU(2) xSU(2) is described in ref. [41] - in that framework, 
the mass of the kaon sets the heavy scale. 

14. As is well known, the straightforward expansion of the masses and 
current matrix elements of the baryon octet |L8| and the vector mesons [42] in 
powers of m u , rrid and m s contains large contributions from infrared singular- 
ities (in particular terms proportional to Mg). It would be of considerable 
interest to apply the reordering of the chiral perturbation series proposed 
in section 0.1 to these quantities: The procedure should lead to a more 
rapidly convergent expansion, so that the first few terms may then provide 
a meaningful determination of the ratio {m d — m u )/{m s — m), for instance. 
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A Low energy representation for the triangle 
graph 

The function j(t) may be expressed in terms of the imaginary part by means 
of the dispersion relation [[[J 

7 (t) = 7 (0) + - / — Im 7 (t') . 



nJ iM 2 t'(t'-ty 

For t — 0, the integral @ can be done explicitly. Expanding the result in 
powers of a = M/m, the subtraction constant becomes 

For a representation of j(t) up to and including first nonleading order, we 
need the imaginary part only to this accuracy. Inserting the representation 
@ and performing the integral, we obtain 

j(t) = ll (t)+ l2 (t) + 0(p) , 

1 f 1 , 2 + JT 2o(2-t) . 2 , 

7i(t) = < — ^= m = H , arcsm 1 — a (ma — 1) 

nyj 32irmM\^F 2-^F tt^t (4-r) 2 

. . 1 fa , / a 
72(*) = ^ ^7 ^=-ln 1 



TT 



32?rmM [ ^4^7 V V^r 

with r = t/M 2 . We have included the subtraction constant in 7i(t), so 
that not only the imaginary part of the term 72 (t), but also its real part is 
significant only in the immediate vicinity of threshold. At leading order of 
the chiral expansion, the expression for 71 (i) reduces to the known result of 
HBXPT ||. The above formulae (i) extend that result to first nonleading 
order and (ii) account for the sum of internal line insertions that determine 
the behaviour in the vicinity of threshold. 

In the notation used in the present paper, the function j(t) represents 
the scalar integral if 21 on the mass shell (q 2 = t, P 2 = m 2 , Pq = 
The above approximate representation for this function is valid throughout 
the low energy region, but is accurate only modulo contributions of 0(p). 
For the calculation of the form factor (section [12|) we need a corresponding 
expression for the infrared part. For that purpose, it suffices to work out 
the regular part i? 2 i to the desired accuracy. As discussed in section || the 
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integral may be represented as (q% = 0, q<i = q, P\ = P) 



R n = ' v " P [ dx / (/: (1 - :) ('- 



r(3 



_- ' / .;,./ .i- i i . \ n 1-3 

(47I")5 Jq Jl 

C = {I- z) {M 2 - x(l - x) q 2 } + z m 2 - z (1 - z) {P - (1 - x) g} 2 . 

The chiral expansions of the regular parts are ordinary Taylor series that 
converge throughout the low energy region. Since the series only starts at 
0(p°), we can truncate it at leading order, i.e. evaluate the above integral 
for M = q = 0, P 2 = m 2 . This gives 

R 21 =m d - 6 2r( d 2 ~^ +0 (p). 

(47T)5(5-d) 

The corresponding counter term was worked out in section [8]. At leading 
order, eq. (|32|) yields z/ 21 = vrC 2 . Indeed the term A v 2 \ removes the pole of 
the T-function at d — 4. Renormalizing at scale fi = m, we obtain 

3 

#21 = -R21 + A z/ 21 = - + 0(p) . 

6l-K l m l 

A corresponding representation for the renormalized infrared part J21 is ob- 
tained by subtracting this term from the formula for ~f(t) given above (for 
an explicit expression, see appendix O). 



B Infrared parts of some loop integrals 

Notation 

m d - 4 [ 1 1 



A 



{^4( hi4 * +r(i > +1 )} 



(4tt) 2 {d-A 
M 



a = — 
m 

1 meson: = iio 

1 f d d k 1 nlf ,/ T 1 , 
A n = - - — -r — — = 2M 2 A + In a 

i Jj (27i) d M 2 -k 2 V !6vr 2 

1 nucleon: An = Ioi 

1 f d d k 1 



A 



N 



i Jj (2vr) a! m 2 - k 2 
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2 mesons: J = 7 2 o 

1 f d d k {1, k», W} 



{J, J", r} = l / 



i A (2vr) d (M 2 - P)(M 2 - (jfe - g) 2 ) 



t = g 2 , 

J{t) = J{t)-2\--^— 2 {2\na + l) 



J( 0s J(0-J(0) = ^|l-^^a«»n^ 

J" = i <f j(0 

JW(*) = {ft - 4M 2 ) J(t) + 2AJ 

J (2 \t) = \j{t)-±- t ^ 

1 meson, 1 nucleon: I = In 

1 /" d d k {1 , A; M } 



(2n) d {M 2 — k 2 ){m 2 — {P — k) 2 ) 
s = p n= s-m 2 -M 2 

2Mm 

JW = i W _£zI2l±«!x, 



r , , i «Vi-n 2 / n + a 

1{ S ) = ~ 7TT i o — Fi 9 arccos 



8tt 2 1 + 2 a n + a 2 \ y/l +2afi + a 2 

1 a(Q + a 



16tt 2 l + 2«H + a 2 



(21na-l) , 



p = p^/(i)( s ) 

/(l)(S) = ^ {( S - ^ + + A 71-} 

2 mesons, 1 nucleon: 

1 A rf d A; {1, k^} 



u in- 1 f d) 



y ( M 2 _ fc 2) ( M 2 _ ( fc _ g)2) ( m 2 _ (p _ A;) 2 ) 

For the low energy analysis of the scalar form factor and of the 7iN scattering 
amplitude to one loop, these integrals are relevant only on the mass shell of 
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the two external nucleons, 

P 2 = P' 2 = m 2 , p> = P - q , t = q 2 . 
The decomposition of the vectorial integral then simplifies to 
I& = {P» + P" t )I& ) {t) + l<fI sl (t), 

4?(t) = 2(4 J_ t) {(2M 2 - t) I 21 (t) - 2/(m 2 ) + 2 J(f)} , 
and the renormalization of the scalar integral reads 

hi{t) = hi{t) H 1 -. arcsm— . 

y/t {Am 2 - t) 2m 

The representation for ^(t) constructed in appendix |S] - which also holds 
near threshold, but neglects terms beyond next-to-leading order - implies 

— , , 1 f 1 , 2 + 1 , / a 



327rmM [ 2 - y/r \ y/I^T 

1 f 2(2-r) v 7 ^ 7T 

H < f arcsin — + ; + 2 lna + 1 ^ + Ofp) 

327r 2 m 2 | ^(4 - r) 2 r 

2 



with r = t/M 2 . 



1 meson, 2 nucleons: 



r 12 



d d k 1 



z J 7 (2vr) rf (M 2 - P) (m 2 - (P - k) 2 ) (m 2 - (P> - k) 2 ) 



In the context of ttN scattering, there are graphs that involve the values of 
I\2 = I\2 {P 2 ? P' 2 ,t) off the mass shell and we therefore leave the variables 
P 2 , P' 2 open. The representation fl28| ) yields 



p 1 /»oo 

J 12 (P 2 ,P /2 ,t) = (4vr)-fr(3-f) / dy dzzC^ 3 , 

Jo Jo 

C = M 2 (l - zf + z 2 {m 2 -y(l-y)t}- 2mMz (1 - z) fi , 

t = (p-p') 2 , n = yn + (i -y)n' , 

n _P 2 -m 2 -M 2 P' 2 -m 2 -M 2 

~ 2mM ' ~~ 2mM ' 

The expression for C shows that the momentum transfer enters exclusively 
through m 2 — y (1 — y) t. The expansion in powers of t thus only yields 
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polynomials - the phenomenon observed in the case of j(t) does not occur 
here. Removing the pole term at d — 4 (see section scaling the variable 
of integration with z-*a« and expanding in powers of a, we obtain 

I 12 (P 2 , P' 2 , t) = I 12 (P 2 , P' 2 , t) + v X2 A , 



v X2 = —^{i- a n- a n' + 0(p 2 )} 

m 



fx(tt) = (1 - a it) Vl -fi 2 arccos (-0) - Q . 

At leading and first nonleading order, the result is independent of the mo- 
mentum transfer. On the mass shell, it reduces to 

I 12 (t) = I 12 (m 2 , m 2 , t) = -^r- 2 (2 In a + 1) + + 0{p 2 ) . 

1 meson, 3 nucleons: 

d d k 1 



'13 : 



% J I (27r) d (M 2 - k 2 )(m 2 - (Pi - k) 2 )(m 2 - (P 2 - k) 2 )(m 2 - (P 3 - A;) 2 ) 



The integral is relevant for ttN scattering. Denoting the momenta of the two 
incoming and outgoing particles by P, q and P', q', respectively, we have 

Pi = P, P 2 = P + q, P 3 = P'-q', P 2 = P' 2 = m 2 , 

q 2 = q ' 2 = M 2 , s = (P + qf = m 2 + M 2 + 2Mm , t = (q- q') 2 . 

The chiral expansion starts at order The coefficients are readily worked 
out with the representation (p8|), setting cf = 4. As in the case of I± 2 , the 
leading and first nonleading terms are independent of t: 



. 1 I 

'13 " 



32vr 2 Mm 3 fi 2 



/o(^) = y/l — Q 2 arccos (-Q) - - - Q . 

C ttN scattering in tree approximation 

In connection with the unitarity condition for the scalar form factor, we 
need the representation of the ttN scattering amplitude to order p 2 . For this 
purpose, the standard decomposition into the invariant amplitudes A and B, 
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is not suitable, because the leading contributions from the two cancel at low 
energies^. We replace A by D = A + u B, 



and work with the variables 



s — u t- 2M 2 



Am Am 
The tree graphs of the Lagrangian + yield 

u y ' > f 2 v 2 -v 2 F 2 F 2 F 2 ^ U \P)> 

B 

B+ (v,t) = 9 -j^^~ 2 +0(p), (C.l) 

„ . , q\m v vu v , o x ^, 5, 

C >>*) = -JTu^-2F 2{9a - 1) + ° {p) 

B 

fi ("»*) = "pT " 2F2 - !) + -pT" + 0(p) ■ 

B 

The projection onto the t-channel I = J = partial wave reads 

o 

Ttl 

/+(*) = + ^^{- 24M2C 1 + ( 4M2 - *) °2 + 6 ( 2M2 - *) C 3}+ °(P 3 ) > 

,om_ ffj^ 3 f arctanr t \ _ - 4M 2 )(4m 2 - t) 



The comparison with the subthreshold expansion of Holder et al. [21 



777/ 

m v B — v 1 2 m z 



10 An evaluation of the scattering amplitude to 0(p n ) suffices to determine D(v, t) to the 
same order and yields B{v, t) to 0(p n ~ 2 ), but fixes A{v, t) only to 0{p n ~ 1 ). Conversely, a 
representation for the pair (D,B) to accuracy (p n ,p n ~ 2 ) fully characterizes the scattering 
amplitude to 0(p n ), while an analogous representation for the pair (A, B) does not suffice, 
even if A is assumed known to all orders. 
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implies the following representations^] for the coupling constants |38 



Cl 

c 3 



F " {d+ + 2M 2 d+} + O{M), c 2 = ^-dt + O(M) 



AM 2 
-F 2 d, 



(C.3) 



The numerical values d. 



1 



C4= 2^ {F26oo "^ + 0(M) - 



oo 



-1.46, df = 1.12, t#i = 1.14, 6, 



7 oo 



10.36 (in 



pion mass units, pi]) thus lead to 



c\' = —O.QOmj/ , <4 = 1-6 , Cg " 1 = —3. 4m JV 1 , C4 = 2.0 m^ 1 . 

The current algebra formula a = -4ciM 2 then yields a = 50MeV. Al- 
ternatively, replacing the input for dJo by d = 45MeV [20], we obtain 



jo) 



— 0. 54 m^ 1 and dtc 



00 



-1.54. 



D Contributions generated by the A 

Effective Lagrangian 

The A degrees of freedom may be described in terms of a Rarita-Schwinger 
spinor ip^, which transforms with the representation x D 1 of the isospin 
group. The condition T a ip® = eliminates the component with I — To 
lowest order, the relevant effective Lagrangian reads 

C A = ^A^C + 9aF < V + h.c.) , (D.l) 

with u a ^ = I tr {r a u^) . The kinetic term involves the tensor 

A"" = -{ip- m A )g^ + i(^D v + j w D") - j"{ip+ m A )Y , 

while the one occurring in the interaction is of the form 

= g, v - (| + Z) lillv . 

We use the notation of ref . |21| , where the coupling constant g A carries the 
dimension of an inverse massQ. The expression for the width of the A that 
follows from the tree graphs of the above Lagrangian reads 

r A -wv. = {(m A + m) 2 -M 2 } , 

247T m A k J 

11 For a discussion of the corrections of order M, we refer to section p"2|.5. 
12 In the literature, the coupling constant <?a is often replaced by the dimensionless 
quantity F„g A = C/V2 §§ |§ = h A § = 97rNA @- 
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where q& is the momentum of the decay products for the decay at rest. In 
the numerical evaluation, we use m& = 1.232 GeV and the value of the cou- 
pling constant advocated in ref. J2I] g& = 13.0 GeV -1 , F n g& = 1.2. There is 



considerable scattering in the numerical values used in the literature. In par- 
ticular, the experimental width implies stronger coupling - the significance 
of the difference is discussed in detail in ref. on the basis of an analysis 
of the P33-wave. Some of the values found within the small scale expansion 
instead indicate a weaker coupling (see ref. |5T| for a review of that expansion 



in HBXPT and ref. || for the corresponding relativistic formulation). 
Contribution to the ttN scattering amplitude 

In addition to the coupling constant g& that determines the residue of the 
poles arising from A-exchange, the above effective Lagrangian contains the 
parameter Z, which generates polynomial contributions to the scattering am- 
plitude (the representation obtained from the tree graphs of this Lagrangian 



can be found in ref. pl| ). To fix the value of Z, we observe that the unitarity 
bounds |$9| ensure an unsubtracted dispersion relation for B + (u, 0)/V. The 
relation represents this amplitude as a superposition of contributions aris- 
ing from the various intermediate states. While the nucleon generates the 
Born term in eq. ( p.l| ), the A produces two peaks in the imaginary part, at 
v = ±co>a, where — { m \ — m2 ~ M 2 )/2m. In the narrow width approx- 
imation, the corresponding contribution to the dispersion integral is of the 
form 



The tree graph amplitude generated by the above effective Lagrangian does 
contain such pole terms, but involves an additive constant proportional to 
Z 2 . In the language of the dispersion relation, this term corresponds to the 
integral over those intermediate states that remain when the nucleon and the 
A are removed. In order for the effective Lagrangian in eq. ( p . 1|) to properly 
account for the contributions generated by the A, we must set Z = 0. 



The corresponding expressions for the invariant amplitudes read [21 



18m L V -1 / V^A _ v + v J \P3 

with z/a = l<ja + t/4m. The coefficients are fixed by the masses. The expres- 
sions may be simplified by using the parameter = (m 2 A +m 2 — M 2 )/2m& 
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(energy of the nucleon in the rest frame of the A): 

«i = 2 (E& + m) {Ea (2 mA + m) — m (r/iA + 2 m)} , 
a 2 = 3 (m A + m)/2 , 

a 3 = —4 m (mi + m) (2 m A + m m A — m + 2 M )/ m A , 
/3i = 2(£ A + m)(£ A -2m) , 
ft = 3/2 , 

/?3 = 2 m(mA + m) 2 jm\ . 

For the reason given in appendix O, the amplitudes A 1 * 1 are not suitable 
to sort out the low energy structure. The two terms in D = + u B ± 
nearly cancel also in the present case: The chiral expansion of D£ and 
only starts at 0{p 2 ) and 0(p 3 ), respectively We do not list the corresponding 
explicit representations - these are readily obtained from the above formulae. 

Subtraction constants and saturation 

The QCD Lagrangian does not contain parameters that explicitly refer to 
the A, nor is this state singled out by chiral symmetry. We are identifying 
the A-contributions on the basis of a saturation hypothesis: The argument 
used to pin down these contributions relies on the fact that the dispersion 
relation for B + {v ) t)/v does not contain a subtraction - we have determined 
the contributions from the A by saturating this relation for t = 0. We now 
briefly discuss the dispersion relations obeyed by the three other amplitudes 
from the same point of view. More specifically, we wish to show that 

• The dispersion relation for D + (u,t) contains a subtraction, which en- 
sures that this amplitude has an Adler zero in the vicinity of the point 
v = 0,t = 2M 2 . 

• The amplitude D~(v,t)/v obeys an unsubtracted dispersion relation, 
which is closely related to the Adler- Weisberger relation and allows us 
to identify the A-contribution to the coupling constant g^. 

• Finally, A~(v, t)/v also obeys an unsubtracted dispersion relation. The 
representation for this amplitude obtained by saturating the dispersion 
integral with the A is identical with the one given above. This confirms 
that the parameter Z must be set equal to zero. 

For simplicity, we only consider the forward dispersion relations, i.e. set t = 
(in the general case, the subtraction "constants" are functions of t). As the 
unitarity bounds of ref . are too weak to discuss the subtractions occurring 
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in the dispersion relations for the isospin odd amplitudes, we invoke Regge 
pole analysis, which yields more specific information about the asymptotic 
behaviour. 

In the case of the isospin even amplitudes, the asymptotics is dominated 
by Pomeron exchange: A + (v,t), D + (v,t) oc v ap< t\ B + (v,t) oc i/^C*) -1 , with 
a>p(0) = 1 (we disregard the logarithmic terms arising from the fact that the 
Pomeron does not represent a simple Regge pole). Since the Regge behaviour 
is consistent with the unitarity bounds, it merely confirms that B + (u, 0)/V 
obeys an unsubtracted dispersion relation. The dispersive representation for 
D + (i>, 0), however, requires a subtraction (note that D + (u,t) is even in u, so 
that a single subtraction suffices). In fact, if this amplitude were to obey an 
unsubtracted dispersion relation, the saturation with the A would lead to a 
conflict with chiral symmetry, which implies that D + contains an Adler zero. 
In the above explicit representation, the zero arises because the subtraction 
constant compensates the pole terms: D^(0,2M 2 ) represents a term of 
order M 4 . The occurrence of a subtraction is essential for this to happen. 

For the isospin odd amplitudes, the p-trajectory yields the leading term: 
A~(v,t),D-(v,t) oc is ap{t) , B~(v,t) oc z/*^*)- 1 , with a p (0) ~ \. This implies 
that the functions 

Mv) = - A~(u, 0) , D(u) = - D~(v, 0) 

V V 

obey unsubtracted dispersion relations (the same also applies for B~(v, 0), 
which represents the difference between the two). Indeed, the representa- 
tion given for A^(y, 0) is precisely the one obtained by saturating the un- 
subtracted dispersion relation for A{y) with the A. The one for D^(u, 0), 
however, does contain a subtraction term, (3 3 , in apparent conflict with the 
statement that this amplitude obeys an unsubtracted dispersion relation. 



Adler- Weisberger relation 

To resolve the paradox, we consider the dispersion relation for D(v), which 
according to the above is free of subtractions: 

A/ \ g\m vb_ 2 f°° dv' , 



F 2 v\ — V 2 TC J M v' 2 — V 2 



We have explicitly indicated the contribution from the one nucleon inter- 
mediate state, with z/g = —M 2 /2m. The comparison with the low energy 
theorem ( |(J.1D leads to an on-shell variant of the Adler- Weisberger relation: 



2F 2 



Od"l) + - r^lmD-(u',0) = O(p"). 

71 JM v 
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The imaginary part may be expressed in terms of total cross sections. Drop- 
ping the corrections^ of order p 2 and replacing F by F n , the result obtained 
for the dispersion integral yields a measurement of g^, based on total ttN 
cross sections. The result is in remarkably good agreement with the value 
found from neutron decay . 

Replacing the dispersion integral by the contribution from the A, the 
above formula yields 

2 2 2 2 {(m A + mf - M 2 } 2 {(m A - m) 2 - M 2 } 

9 A + g * 9A ^ ^ _ m 2 _ M 2)2 

The corresponding numerical value, qa = 1.35, shows that the contribution 
from the A indeed dominates the dispersion integral: The remainder is about 
four times smaller (and of opposite sign). 

This resolves the paradox encountered above: The explicit representa- 
tions for the amplitudes B A do not include the effect generated by the 
A in the value of the coupling constant g^, which is treated as an indepen- 
dent parameter - this is why the representations for D A and B A involve the 
subtraction term /3 3 , despite the fact that the full amplitudes obey unsub- 
tracted dispersion relations. The difference only concerns these two ampli- 
tudes: The formulae (|C.1| ) show that the polynomial term proportional to 
g\ — l, which causes the paradox, only enters the expressions for D~ and B~ 
(in A~ = D~ — v B~ , the term drops out). 



Subthreshold expansion 

The A-contributions to the coefficients of the subthreshold expansion are 
obtained by expanding the amplitudes in powers of v and t, with the result 

io moj/\ y muj A 

rf 01A = ~ i Q 9 \ 2 ("1 - 4 "2^ A ) , &00A = ~TF^ ( 2 A ~ U M ■ 

Numerically, this amounts to (c^ooA; ^ioa> ^oia> ^ooa) = ( — 1-4, 0.8, 0.7, 5.3), 
in units of M n . The comparison with the phenomeno logical values of ref. |2T| , 
(c/q , di Q , doi) &oo) = (~ 1-46, 1.12, 1.14, 10.36), shows that the A indeed ac- 
counts for a significant part of the curvature seen in the scattering amplitude 
at small values of v and t. The situation is similar to the one in the mesonic 
sector, where the most important singularity that is not explicitly accounted 
for in the effective theory is the one generated by the p. 



13 



For an analysis of the terms of 0(p 2 ), we refer to Brown, Pardee and Peccei Eol 
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It is instructive to compare the chiral expansion of the subthreshold co- 
efficients with the small scale expansion thereof. In the chiral expansion, 
the coefficients are expanded in powers of M at fixed m^. The small scale 
expansion results if the mass difference A = tua — m is also treated as a 
small quantity and the double series in M and A is reordered by counting 
the ratio M/A as a quantity of order 1. The chiral expansion yields a decent 
approximation for all of the above coefficients (the largest deviation occurs in 
df Q , but it still amounts to less than 20 %). The small scale expansion does 
not improve the accuracy of this approximation - in fact, the discrepancy 
between the leading term of that expansion and the full result for df is about 
twice as large as for the chiral series. The numerical exercise shows that - 
as long as we stick to small values of v and t - the singularities generated by 
the A do not jeopardize the chiral perturbation series. 

Contribution to the effective coupling constants of Cff 

To read off the contributions of the A to the effective coupling constants of 

(2) 

C N , we expand the above expressions in powers of M and insert the result 
in eqs. (|3"9"|). This leads tdQ. 

cf = , = 2cm (m A + m) , cf = — c (4m A — wa™ + m 2 ) , 

a 2 F 2 

cf = c(2m& - m)(m A + m) , c = - — ^—^ -. (D.2) 

9m A (m A — m) 

Numerically, these formulae yield (cf , cf , cf , cf) = (0, 1.4, —2.0, 1.1) m^ . 
Again, the higher order terms of the chiral expansion do not significantly 
modify these estimates: Evaluating the relations ( |3"9"D with the full expres- 
sions for the coefficients of the subthreshold expansion instead of only the 
leading terms, we obtain (cf , cf , , cf) = (—0.04, 1.2, —2.1, 1.2) m^ 1 . In 
particular, the higher order effects generate a small contribution to the cou- 
pling constant c±, which amounts to a shift in the a term by about 3 MeV. 



Contribution to A a and to the coupling constants ei, e2 

As discussed in section 0.4, the corrections to the leading order result for 
the difference A a = a(2M 2 ) — er(0) are dominated by the logarithmic term, 
A2 M 4 InM/m. A significant fraction of this term is generated by the A: 
Inserting the above expressions in eq. (|38|) , we obtain 

g 2 A (Qm 2 A -2m A m + m 2 ) 
l A 367T 2 m A (m A — m) 



14 Note that the coupling constant C4 picks up a small negative extra contribution from 
the term c£ = — l/4m in eq. (|C.3|). 
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Numerically, this amounts to A2 M 4 In M/m = 5.9 MeV, to be compared 
with the number 10.6 MeV, obtained by inserting the phenomenological val- 
ues for the combination C2 + 6 C3 of effective coupling constants (the contri- 
bution from gA is of opposite sign and reduces this to 7.8 MeV). 

The effects generated by the A in the scalar form factor were analyzed 
quite some time ago p9j . That calculation is equivalent to an evaluation 



of the dispersion integral flUP in the approximation where the integrand is 
replaced by the leading term of the small scale expansion. The corresponding 
representation for the contribution to the partial wave /?(£) generated by A 
exchange reads 



p0M 2g\m 2 (2A 2 + t-2M 2 \/t — A M 2 

TA4.it) = < : — arctan z 

M+W 9tt \ Vt-AM 2 2A 

with A = rriA. — m. At low energies, this approximation yields a remarkably 
good representation for the full partial wave projection of the A-contribution 
to the scattering amplitude. The leading term in the chiral expansion of the 



dispersion integral then reads [|32 



g\M' [5 tt (M\ Mn fM 2 

The coefficient of the chiral logarithm occurring here indeed agrees with the 
leading term in the small scale expansion of the representation (p.3| ). Note 
that there is a difference between the chiral perturbation series and the small 
scale expansion: In the former, chiral logarithms only start showing up at 
0(p 4 ), but for the latter, those generated by the A manifest themselves 
already at third order. 

The remainder of the above representation amounts to an estimate for 
the A-contribution to the coupling constant e 2 : 

A g'i f , /2A\ /M 2 
~ A A J 10 1n +7 + — 



2 144tt 2 A { \m J \A 

The formula confirms that the coupling constant e 2 is small: Numerically, 
the corresponding contribution to A CT amounts to 0.7 MeV. 

The analogous estimate for the A-contribution to the low energy constant 
e\ may be obtained from the small scale expansion of the nucleon mass |32 
The comparison with eq. (0) yields 



48tt 2 A [ \2mJ \A 2 
Numerically, this corresponds to a contribution to the a-term of 5.7 MeV. 
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Note that, if the degrees of freedom of the A occur as dynamical vari- 
ables of the effective theory, the quantum fluctuations also generate graphs 
that contain A propagators instead of nucleon propagators. The small scale 
expansion allows a coherent counting of powers also for these graphs. The 
width of the A, for instance, represents a term of third order, so that the 
corresponding shift of the pole in the propagator does represent a correction 
which perturbation theory is able to cope with. There is a difference to the 
framework described above, insofar as the constants m and c\ now pick up 
renormalization. At leading order, the A is degenerate with the nucleon, 
so that the loop graphs containing A propagators also give rise to infrared 
singularities. The method described in the present paper may be extended 
to analyze these ||. 
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